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The Effect of Periodic Pressure and 
Temperature Fluctuations on Unsteady 
Heat Transfer in a Closed System 

DANIEL W. WENDLAND 
SUMMARY 

The purpose of this project was to study the nature of the 
unsteady heat transfer from a closed gas system to a plane 
metal surface with the gas experiencing periodic pressure 
and temperature oscillations. The gas system was defined 
by a pis ton- cylinder-head compression chamber; the piston 
was driven by a variable-speed motor to generate the gas 
pressure and temperature fluctuations. Gas pressure his- 
tories were measured with a piezoelectric pressure trans- 
ducer. Gas temperatures were measured with a 3.81 micron 
hot-wire anemometer probe used as an uncompensated resis- 
tance thermometer. Instantaneous metal surface (the sur- 
face of the head of the chamber) temperature was measured 
with a fast-response surface thermocouple with a one-micron 
junction thickness. 

A mathematical technique was developed to calculate in- 
stantaneous metal surface heat transfer rates: experimen- 

tally measured surface temperature histories were intro- 
duced into a CDC 3600 digital computer program which cal- 
culated, listed, and plotted surface heat transfer his- 
tories . 
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In the experimental observations, the gas pressure history 
had an amplitude of about 365 psi and a mean value of about. 
197 psi; the pressure maximized at about 2 1/2 degrees be- 
fore the volume minimum (TDC). These histories were essen- 
tially unaffected by system frequency throughout the exper- 
imental range (15 to 58 Hz). The bulk gas temperature his- 
tories had amplitudes on the order of 600 P° , mean values 
near 370°F, and peaked at about three degrees before TDC; 
these histories were somewhat frequency-sensitive. Some of 
the frequency-sensitivity is attributed to phase lags of 
the gas temperature transducer. 

As the driving frequency was changed, two regimes of re- 
sponse in periodic metal surface temperature to these gas 
property variations were noted. At moderate and high fre- 
quencies (25 to 58 Hz) the wall surface temperature was 
fairly sharply peaked, maximizing at about zero crankangle 
degrees (i.e., piston TDC). A theoretical model viewing 
the thermal boundary layer as a time-varying resistance 
predicted essentially this behavior. An even better pre- 
diction was obtained with a model derived from the one- 
dimensional energy equations and containing resistive, ca- 
pacitive, convective, and pressure work terms. These models 
are not included in this report, as a far superior theoret- 
ical model was developed in the last months of the research 
period. For information on these secondary models, the 
reader is referred to the appendices of the author’s Ph.D. 
thesis, published in January 1 , 1968 at the University of 
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Wisconsin. 

The superior theoretical model mentioned above was formu- 
lated from thermodynamic principles and contained little 
or no empirical information. The model predicted wall 
heat transfer histories which were very similar to the heat 
transfer histories calculated from the medium and high fre- 
quency experimental data. The model, termed the Adiabatic 
Plane Model, also predicted negative, zero, and infinite 
heat transfer coefficients, and a phase lag of wall surface 
temperature with respect to bulk gas temperature of less 
than one crankangle degree. 

At very low frequencies (up to perhaps 18 Hz), the wall 
temperature response was smooth and nearly sinusoidal, lag- 
ging the gas temperature peak by as much as 80 crankangle 
degrees; the peak location was highly frequency-sensitive. 

A theoretical model for the boundary layer proposing only 
resistive and capacitive effects can predict this type of 
temperature response. This model is explained in the 
appendices of the thesis mentioned above. 

At intermediate frequencies, the two regimes combine and 
the wall surface temperature shows two maxima per cycle. 

In summary, three conclusions can be drawn from this in- 
vestigation: 

1) the metal surface temperature and surface heat 
transfer rate are not in phase with the gas tern- 


perature because of compression and expansion of 
the gas; the heat transfer coefficient based on 
Newton's Law of Cooling is therefore of doubtful 
utility, as the coefficient assumes negative, 
zero, and infinite values; 

2) the lag of surface temperature behind gas tem- 
perature and pressure histories at moderate and 
high frequencies was found to be small (on the 
order of two or three crankangle degrees , or 
about .04 radians) and insensitive to frequency; 

3 ) the surface heat transfer rate was found to lead 
the gas temperature at moderate and high fre- 
quencies by from eight to ten crankangle degrees 
(0.18 radians). 
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INTRODUCTION 


One of the classical problems in the study of heat trans- 
fer is the case of energy exchange between a fluid and the 
plane surface of a solid plate or wall. Typically, the 
fluid at a great distance from the plate is at a constant 
temperature, and the plate is at a different temperature. 
The fluid traditionally moves in steady flow parallel to 
the plate, heat transfer occurs between plate and fluid, - 
and velocity and temperature profiles -are established. - 
in the fluid. In laminar flow, conditions at a given 
point of reference with respect to the solid surface do 
not change with time. In turbulent flow, relatively small 
perturbations in fluid properties are superimposed upon 
the unchanging flow conditions . Steady-state heat trans- 
fer at the plate surface at a distance x from the leading 
edge of the plate is described by: 

q(x)/A - h(x) • [T g - T w ] (1) 

where: h(x) is the local heat transfer coefficient, 

T g is the local plate surface temperature, and 
T*, is the fluid temperature at a great distance. 

Theoretical arguments and experimental studies have pro- 
duced many correlations for h in terms of fluid properties 
and the fluid velocity (or Nusselt Number in terms of 
Prandtl Number and Reynolds Number) . 



This investigation, too, is concerned with energy trans- 
port between a fluid and a wall. There are several con- 
ditions in this study, however, which differ significantly 
from the "classical" conditions: 

1) The fluid motion is normal to the wall, not 
parallel to it. 

2) The system is closed, not open; one mass of 
fluid (in this case, gas) undergoes a series 
of thermodynamic processes. 

3) The gas conditions (pressure, temperature, 
velocity) far from the solid surface (in fact, 
throughout the system) are cyclical with time. 

4) The plate, or wall, conditions (surface temper- 
ature, surface heat flux) are also cyclical. 

This type of gas-solid interaction is of interest ir. two 
general areas: the study of unsteady heat transfer in 

rocket combustion chambers, and research in internal com- 
bustion engines (and in compressors). 

In the piston engine, much of the energy of combustion and 
compression is lost through heat transfer. The precise 
nature of this unsteady heat transfer is still largely un- 
known; steady-state measurements (i.e., cooling water tem- 
peratures and flow rates) tell nothing about the time dis- 
tribution of the energy loss from the gas around the work- 
ing cycle. In a motored engine, the gas temperature may 
change by six or seven hundred Fahrenheit degrees during 
the cycle, while the response in compression chamber wall 
surface temperature will be on the order of a few degrees. 
Temperature gradients in the gas immediate to the wall may 


be as high as 15,000 P°/inch, the time rate of change of 
gas temperature may be as high as 2.5 x 10^ F°/second 
(12 F° per crankangle degree at 3600 rpm), and the time 
rate of change of gas pressure may reach 1.3 x 10^ psi/ 
second. 

In the rocket combustion chamber and nozzle, an oscilla- 
tory combustion phenomenon often occurs in which the gas 
pressure can swing by several hundred psi and the gas 
temperature by over 1000 F° at frequencies from 50 - 250 
Hz ("chugging") to several kHz ("screaming"). This rep- 
resents rates of change of gas temperature and pressure 
several orders of magnitude greater than those encountered 
in engine work. Heat transfer takes place from hot gas to 
the walls; when the gas - wall temperature gradients are 
high and the thermal resistance of the gas near the walls 
is low, the heat flux rates are often high enough to burn 
through the nozzle. 

Heat transfer also takes place from the hot gases to the 
vaporizing fuel and oxidizer droplets near the injector. 
The vaporization rate of the droplets responds in some 
fashion to the gas pressure and temperature changes; an 
increase in vaporization rate may support or oppose the 
impinging pressure waves, depending on the phase rela- 
tionships . 


The thermal response of surfaces to changes in gas con- 


ditions is complicated by two factors: 

1) as the g. has a heat capacity, capacitive energy 
storage ana release will take place in the gas 
immediate to the wall at a rate proportional to 
the time rate of change of local temperature, and 

2) since the gas pressure is changing with time, 
work done by the volume of gas near the wall will 
represent a source or a sink of energy. 

Thus the behavior of the temperature and the energy fluxes 
in the wall and in the gas near the wall surface is of 
great interest. This is a broad area for inquiry: the 

cyclic frequencies can range .to 10,000 Hz., the energy 
fluxes to hundreds of thousands of Btu/hrft 2 , the temper- 
atures to several thousand degrees R. The medium can 
range from air to luminous combustion gases. It is also 
a very difficult area for inquiry: rapidly changing local 

temperatures must be measured in complex flow situations 
under conditions far from ideal. 

Practical considerations make it necessary to limit an 
investigation to a small range in operating variables and 
flow conditions. This investigation limits itself to fair- 
ly low frequencies (as far as rocket combustion instability 
is concerned), moderate pressure and temperature ratios 
(or swings), and a closed gas system. A study concerned 
with higher frequencies (= 500 Hz), smaller pressure and 
temperature ratios, and an open system will be published 
in 1968 by R. W. Goluba*. 

*Ph. D. candidate. University of Wisconsin ( 1968 ) 
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I. SURVEY ON PRESSURE WAVES AND HEAT TRANSFER 

Prior to 1950, little interest was shown in the influence 
of pressure waves, mechanical vibration, and fluid velo- 
city pulsations on heat transfer. Since that time, how- 
ever, and particularly in the last ten years, quite a bit 
of effort has been expended in this direction. This is 
probably due to a number of factors: the discovery that 

certain frequencies and amplitudes of mechanical vibration 
can improve the performance of heat exchangers, the occur- 
rence of vibration in rocket combustion and heat transfer, 
and the recent availability of high frequency-response re- 
cording equipment to study these effects. 

I V . . 

Interest has centered, largely, on two types of systems. 

t , 

In the first, a horizontal or vertical plane surface or 
small cylinder in a fluid environment is mechanically in- 
duced to vibrate, and the influence of this motion on free 
or forced convection from object to fluid is noted. Typi- 
cally, the heat transfer rate changes little until a criti- 
cal vibration intensity (the product of frequency and am- 
plitude) is reached, at which point the heat transfer rate 
increases by some 50 %. The critical vibration intensity 
generally correlates with a transition in the type of boun- 
dary layer adjacent to the surface. Shine (1) presents a 
good example of this work and a bibliography of the same. 


The second, and more pertinent, system category is one in 
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which pressure and/or velocity pulsations in a fluid en- 
vironment impinge on a stationary surface (horizontal or 
vertical, plane or cylindrical), and the rate of heat 
transfer in free or forced convection is measured as a 
function of pulsation frequency and intensity. These re- 
searches can generally be categorized by the intensity, or 
amplitude, of the disturbance. 

The great preponderance of papers found in the literature 
deal with very low amplitude (under one psi), or acoustic, 
oscillations in pressure and the attendant pulsations in 
fluid velocity. A very few investigations deal with high 
amplitude (up to, say, 100 or 200 psi) disturbances which 
are generated by sirens, resonators, or mechanical devices. 

A small but increasing number of studies are being made 
with very high amplitude (on the order of 500 psi) pressure 
waves whose source is an oscillatory combustion. 

The low amplitude, or acoustic, disturbances have received 
much attention. Little theoretical work has been done in 
the area; most investigations are directed toward the mech- 
anisms of heat transfer in the boundary layer with and 
without acoustic stimulation. Early investigators were 
concerned primarily with flow patterns in resonant tubes 
and not with heat transfer. Lord Rayleigh (2) performed a 
theoretical analysis of the flow patterns in Kundt's tube 
and predicted that resonant vibrations would produce quarter 
wave length vortices. This was later experimentally veri- 
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fled by Andrade (3). C. K. Hwu (4) investigated the effect 
of acoustic vibrations on the forced convective heat trans- 
fer for air flowing in a horizontal tube. An increase of 
up to 50 % in laminar flow heat transfer was noted when the 
tube length and acoustic frequency were such that a stand- 
ing wave was set up in the tube. 

Pand (5) has done much work with a heated horizontal cylin- 
der in an intense sound field. Free convection from the 
cylinder was found to be invariant with sound intensity 
below about 135 db (i.e., peak to peak pressure swing of 
0.046 psi); from 136 db to 142 db at 1 kHz to 6 kHz driving 
frequency (i.e., 0.052 psi to 0.103 psi peak to peak), the 
rate of heat transfer increased sharply. Above 146 db 
(0.164 psi peak to peak), the increase was less rapid. 
Harrison, et . al., (6) studied the influence of stationary 
sound waves on a vertical heated tube in both free and 
forced convection. Up to 542 increase in film coefficient 
was observed. Kubanskii (7,8) has presented a number of 
excellent papers in this area. An electrically heated hor- 
izontal cylinder was subjected to an intense sound field 

i 

during free convection in air. With frequencies from 8 kHz 
to 30 kHz, a 75 % increase in heat transfer rate was noted 
for the intensity range of 145 db to 152 db (0.146 psi to 
0.326 psi peak to peak). The same horizontal tube in forced 
convection was exposed to standing sound waves at right 
angles to the tube. The heat transfer was either not 
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changed or increased by about 50% t depending on the re- 
lation between the nodes of the standing wave and the tube 
axis. A practical application of this increase in heat 
transfer was noted by Boucher (9), who reported that the 
application of sound waves of 1^5 db intensity and 6 kHz 
frequency greatly increased the drying rate of liquid 
sprays, raw cotton, rubber pellets, and many other materials. 

Researchers at the Georgia Institute of Technology have 
done much work in the area of acoustics and heat transfer 
since 1958 under the sponsorship of the Aeronautical Re- 
search Laboratory at Wright Field. Their work is summar- 
ized in papers by Jackson (10), Purdy (11), and Eastwood 
(12). The outside surface temperature of a 10.5 foot, 

^•125 inch outside diameter horizontal copper tube is 
maintained at 210°F by a jacket of saturated steam. Air 
flows through the tube at about atmospheric pressure. A 
resonant acoustic vibration is generated at the downstream 
end by a 50-watt driver (or horn) operating in the fre- 
quency range 150-350 Hz. The greatest sound intensity ob- 
tained was 162.5 db (approximately 1.1 psi peak to peak) at 
222 Hz. 

Experimental results indicate that the resonant sound waves 
produce a spatially periodic variation in the local heat 
transfer coefficient. Two distinct spatial distributions 
of the heat transfer coefficient were discovered as the 
Reynolds number was changed. For Reynolds numbers below 
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about 20,000 the local heat transfer coefficient maxi- 
mized at the velocity antinodes (which are also the pres- 
sure nodes). Above a Reynolds number of about 40,000 the 
heat transfer maxima are at the velocity nodes (which are 
the pressure antinodes). The intermediate Reynolds numbers 
produce a gradual transition between the two types of heat 
transfer, the actual changeover taking place at Re=35,000. 

A locus of critical sound intensity levels exists below 
which essentially no change in the heat transfer takes 
place. The critical sound intensity level also shows two 
distinct zones of behavior with Reynolds number; in each 
zone, the critical intensity level is directly proportional 
to Reynolds number. For the low Reynolds number case, an 
increase in sound intensity level above the critical value 
produces a slight (at most, 20 %) increase in the overall 
heat transfer level. For the high Reynolds number case, 
increases in intensity above the critical tend to decrease 
the overall heat transfer level. The reason for the shift 
in behavior with increasing Reynolds number is ascribed 
to changes in secondary flow patterns. 

Other papers describing experimental and analytical inves- 
tigations of the influence of low— amplitude pressure waves 
on heat transfer are given by Westervelt (13)’, Illingworth 
(14), and Lemlich (15), to list a few. 

Large amplitude (50 psi to 150 or 200 psi peak to peak) 
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steep-fronted pressure waves are being studied by several 
groups. At the University of Wisconsin, R. W. Goluba (1 6) 
is working with siren-generated pressure waves in a one- 
inch diameter tube. The upstream siren is a rotating disk 
which chops the flow and produces a resonating, steep- 
fronted pressure wave at 500 to 1100 Hz. The mean tube 
pressure is from 200 psi to 300 psi; the peak to peak pres- 
sure amplitudes generated are from 100 psi to 120 psi. 
Transient tube wall surface temperatures are very similar 
to the gas pressure histories in shape, from 0.25 to 0.75 F° 
in amplitude, and generally in phase with the gas pressure. 
Generalized correlations for heat transfer have not yet 
been obtained by Goluba. 

At Princeton University, studies under the direction of 
Crocco (17) of large-amplitude waves have been in process 
for nearly ten years. The oscillations (about 50 psi peak 
to peak amplitude) are steep-fronted: the disturbance 

caused by a rotating chopper located downstream of the test 
section has steepened into a shock front by the time it 
reaches the test area. Heat transfer increases of up to 
200# have been observed, and are generally associated with 
flow reversals; recall that secondary flows also played an 
important role in the Georgia Institute of Technology in- 
vestigations. Although no generalized correlations have 
been developed by the Princeton group, the increase of heat 
transfer with resonant pressure oscillations has been found 
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to vary with oscillation frequency and amplitude, the 
through-flow Mach number, and the mean tube pressure. The 
local oscillating pressure, however, was shown not to be 
the driving parameter for the heat transfer increase; the 
amplitude of the local oscillating velocity provided a much 
better correlation. These large amplitude oscillations are 
normally associated with flow reversal. The Princeton 
group has had difficulty separating amplitude effects from 
secondary effects due to wave shape changes and deforma- 
tions in radial velocity profiles. 

Many investigations, both analytical and experimental, are 
currently being conducted to determine the driving mechan- 
isms and suppression techniques for large amplitude, steep- 
fronted pressure waves in oscillatory combustion. Much of 
this information remains in the classified literature. 

Among the experimental contributions, Blackshear (18) 
presents a good, basic study of "screech'' mechanisms in a 
small afterburner. Barrere (19) has published a lengthy 
experimental examination of low frequency combustion in- 
stability including reproductions of instantaneous gas 
pressure and temperature recordings. Cheng (20) of the 
Guggenheim Laboratories authored an experimental study of 
small amplitude, low frequency pressure oscillations in the 
combustion chamber of a liquid-fueled rocket with different 
nozzle profiles. Marble (21) traces "chugging" in a bi- 
propellant system to fuel and oxidizer flow rate oscilla- 



tions and presents a servo-stabilization mechanism to 
sense and reduce the "hunting” of the system. 
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Schlieren photographs of high frequency oscillations in a 
bluff-body flame are given by Kaskan (22) and are associated 
with a transverse oscillation mode. Truman, et . al., (23) 
has also done a study of the transverse mode of vibration, 
but in high thrust engines. Finally, Rogers (2 1 *) and Ber- 
man (25) have published experimental studies of high fre- 
quency instability, the former in ramjet combustors and 
afterburners, the latter in rocket nozzles. The publica- 
tion by Rogers contains good illustrative and background 
material. 

Theoretical studies, primarily on the coupling of pressure 
and velocity disturbances to heat transfer, have been con- 
tributed by a number of authors. Cheng (26), mentioned 
earlier for an experimental study, has presented a theoret- 
ical study on the self-excitation of high frequency combus- 
tion instability in solid propellant rockets. In an earlier 
paper, Cheng (27) elucidates on the effect of non-uniform 
time lags between the introduction of combustibles into a 
system and their conversion to burned gas on the high fre- 
quency instability of the system. 

Crocco, with Cheng, presented a long theoretical disserta- 
tion (28) on the role of unstable pressure oscillations in 
combustion instability with "concentrated combustion." Fi- 
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nally , PriGin (29) has attempted to correlate high frs*” 
quency combustion instability in liquid-fueled rocket com- 
bustors by means of similarity parameters. There are many 
more references to be found in the literature pertaining 
to combustion-related pressure oscillations. 

To conclude this brief review, a number of papers are given 
in which the relationship of flame front energy generation 
to pressure waves is investigated. 

Chu (30) has formulated expressions for the pressure field 
generated by a moderate rate of heat release within a fluid. 
Two general models are considered: in the first, the nature 

of the steep-fronted pressure wave generated by heat release 
in a tube section is determined j in the second, the pressure 
field due to a point energy source in a three-dimensional 
environment is calculated. Blackshear (31) discusses in- 
stability in burners: The dynamics and amplification of 

pressure waves passing through a flame front are analyzed. 

A change in the flame front area is found to act as a gen- 
erator of pressure waves. Schalla (32) corroborates this 
finding: a source of pressure oscillations is present in 

any system in which the gas experiences an increase in heat 
release rate.. Using ethylene-air mixtures, Schalla found 
acoustic pressure disturbances (on the order of 0.10 psi 
amplitude, 120 Hz frequency) generated when the air flow 
was across sharp gradients in composition. 


The amount of publication in the area of heat transfer 
augmentation due to periodic pressure and velocity oscilla- 
tions can be expected to increase in the future. 


II. EXPERIMENTAL APPARATUS 
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GENERAL STATEMENT ON APPARATUS 

The object of this investigation is to examine the effects 
of periodic pressure and temperature oscillations in a 
closed gas system on the heat transfer between the gas and 
a metal wall. A piston-cylinder-compression chamber com- 
bination is used to generate the periodic fluctuations in 
pressure and temperature. This compression chamber was 
taken as a reasonable source of low-frequency fluctuations 
for a number of reasons : 

1) A pressure ratio as high as 25 could be produced, 

\ and it could be easily changed by changing the 

compression ratio of the system; 

2) A broad low-frequency range of 9 to 60 cps could 
be encompassed; 

3) The area available for instrumentation was not 
small; 

*0 The head (i.e., metal wall) material could easily 
be changed . 

The project started with an already set-up CFR engine 
(spark ignition, variable compression ratio), but the in- 
strumentation area was severely limited by the water jacket. 
An air-cooled Briggs and Stratton model 1*12300 engine was 
installed, but oil droplets from the piston rings and cylin- 
der walls coated and/or destroyed the delicate instrumenta- 
tion. A Pioneer model 700 chain saw crankcase was substi- 
tuted: it was designed to need no oil in the crankcase. 
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and the crankshaft and connecting rod bearings could run 
on a dab of grease. To keep the cylinder walls free of 
oil, the steel piston rings were removed, the ring grooves 
were widened and deepened, and teflon piston rings support- 
ed by O-rings were installed. A sprayed teflon film re- 
placed oil as the cylinder wall lubricant. 

The head was a half-inch thick plate of either cast iron 
or aluminum. Wo valves were used, as air turbulence was to 
be minimized; one-dimensional air flow normal to the head 
surface was desired. As blowby occurred past the piston 
rings, makeup air was supplied through two small holes in 
the cylinder wall which were exposed by the piston near the 
bottom of its stroke. 

CRANKCASE AND CYLINDER 

The cylinder bore was 2.3125 inches, and piston stroke 
was 1.5625 inches. Two 1/8 inch diameter holes in the 
honed cast iron sleeve were exposed by the top piston ring 
at 38 crankangle degrees before bottom dead center to supply 
makeup air to compensate for blowby losses. Pig. 1 shows 
the makeup air supply system. The pressure of the makeup 
air system was variable from zero to about 25 inches of 
mercury. 

PISTON 

The piston was cylindrical. Its two main piston rings 
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* 

of teflon were 1/8 inch wide and 0.07^0 inches thick, 
pressed to the cylinder wall by two 0.100 inch diameter 
0-rings. After considerable piston slap was noted from 
the scuffing and dirt build-up on the piston, the piston 
and cylinder were replaced and a thin teflon ring was added 
to the piston skirt. This reduced the scuffing so that 
after eight hours of operation at an average speed of 2000 
rpm, only minor piston scuffing was noted. In nine hours’ 
running, an average of only .003 inches of wear occurred 
on the three teflon rings. 

The cylinder walls and piston sides were sprayed with 

* 

Rocol I. F. L. spray teflon , a new product designed as 
a dry lubricant for die and press work. 

HEAD 

The head was a half- inch thick plate, about four inches 
square, with a one-inch threaded hole in the center for 
temperature measuring instrumentation and a pressure pick- 
up mounting hole slightly off-center. At the cycle fre- 
quencies studied, the head acted as a semi-infinite plate - 
to the temperature fluctuations on the inside surface: 
only steady-state heat transfer occurred at the outside 
surface. The compression chamber side of the head was 
counterbored 0.020 inches concentric with the bore at the 

*See List of Manufacturers and Addresses in rear of text 
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bore diameter to establish a compression ratio of 10.06. 

The off-center location of the pressure pickup necessitated 
the machining of a small, semicircular volume in the head; 
an identical volume was machined symmetrically across the 
bore in the interest of aerodynamic symmetry. These 
volumes (hereafter termed "ears") are visible in Fig. 2. 

The expression used in this work for compression chamber 
volume is: 

-□2 d 2 

V gas * ^clearance + 4 CR(l-coswt) + ^ (l-cos2u>t)] (2) 

" here ^clearance “ °-7253'< inches 3 , 

B = Bore = 2.3125 inches, 

R = Crank Radius « 1.1875 inches, 

L ■ Connecting Rod Length = 4.42 inches, and 
cat = angle referenced to top dead center. 

GENERAL OVERVIEW ON INSTRUMENTATION 

The principal variables measured are gas pressure, bulk 
gas temperature (sometimes simply termed gas temperature), 
wall surface temperature, and driving frequency. Gas 
pressure is measured by a 601H Kistler piezoelectric pres- 
sure transducer. Gas temperature is measured by a small- 
wire resistance thermometer. Wall (i.e., head) surface 
temperature is measured by a surface thermocouple. Driving 
frequency, and crankangle position of the piston (i.e., the 
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time reference), is measured by an Electro* magnetic 
pickup excited by a set of grooves in a crankshaft fly- 
wheel . 

A separate set of breaker points actuated by a cam on the 
crankshaft, with associated circuitry, provides a 10-volt 
square wave at the driving frequency for oscilloscope 
synchronization. The more important of these instruments 
will be covered in succeeding sections. 

ASTRODATA AMPLIFIERS 

While some of the transducers in the experimental setup 
had output voltages on the order of a few volts (timing 
mark pickup, pressure transducer - charge amplifier com- 
bination), the temperature-sensing transducers yielded 
millivolt or microvolt level signals which required am- 
plification before display. Preliminary work in this study 
was done with the hand-built direct-coupled amplifiers de- 
signed and fabricated by Overbye (33); because of stability 
and noise problems, it was deemed necessary to substitute 
a commercial product. Low noise Astrodata* direct-coupled 
amplifiers were obtained. Each amplifier had a gain vari- 
able from 3 to about 3100, differential input and output 
channels, and input and output grounds independent of 
chassis ground. The stated noise figure was less than 2 

•Trade name: see List of Manufacturers and Addresses 
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microvolts rms referred to the input from 0.05 Hz to 
10 kHz. The quoted frequency response was flat ±0.5? from 
dc to 1 kHz, down 3 db at 10 kHz. 

In testing the amplifiers prior to installation, it was 
discovered that they caused a frequency-dependent phase 
shift of the output signal to occur relative to the input 
•signal. The measured mean lag of several amplifiers is 
given in Pig. 3 as a function of frequency with an indica- 
tion of the spread at each measured point. It is apparent 
that below about 5 kHz (i.e., the 83rd harmonic at 3600 
rpm, or 60 Hz), the phase lag for a given amplifier is very 
closely linear. The phase lag could be expressed as: 

. lag of output with 

9 respect to input “ a 

where f is the frequency in Hz and a is the measured slope 
in degrees per kHz. The measured slope over the range 
0 « f $ 5000 Hz varied among the 11 amplifiers tested from 
18.7 degrees/kHz to 19.7 degrees/kHz. Five amplifiers 
with slopes within a percent of 19.3 degrees/kHz were selec- 
ted and used in the project. In all cases the phase shift 
was essentially independent of gain. 

Since the phase lag is linear, the output signal must be a 
faithful (though delayed) reproduction of the input, for if 
the fundamental of the input wave, say 50 Hz, is delayed by 
0.965 degrees, and the seventh harmonic (i.e., 350 Hz) is 


1000 


(3) 



FREQUENCY, KILOHERTZ 




24 

delayed by 6.755 degrees, the time delay of the two is 
identical (0.0536 milliseconds) since the period of the 
higher harmonic is one seventh as long. This does intro- 
duce a complication, however. If a phase-shifted surface 
temperature is displayed on the oscilloscope with a non- 
phase-shifted timing mark trace (since the latter required 
no amplifier), the two traces are obviously displaced hor- 
izontally on the display and cannot be easily related. The 
problem is solved by using Astrodata amplifiers in the 
timing mark and pressure circuits at low gain, and having 
them serve as "phase shifters.". Then all signals displayed 
simultaneously are phase shifted in equal amounts, in- 
cluding the timing-mark base. 

Before the advent of reliable solid state dir.ect coupled 
amplifiers, it would have been very difficult, if not 
impossible, to record or study the small signals studied 
in this project. 




VII . EXPERIMENTAL DATA 
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The experimental data gathered in this investigation con- 
sisted chiefly of temperature, pressure, and heat flux 
histories for the equilibrium cycle of the system using 
air as the working fluid. At each of the eight fixed driv- 
ing frequencies in the range 15 to 58 Hz (915 to 3452 rpm), 
from two to five cycles (not necessarily successive cy- 
cles) of data were taken; each cycle of data consisted 
primarily of the bulk gas temperature and pressure his- 
tories and the head surface temperature history. The in- 
strumentation used for these measurements has already been 
described. A table listing the figure numbers for most of 
the experimental measurements is given below: these figures 

are found in Appendix A; data for three intermediate speeds 
are not included. 

Pigure Numbers for Experimental Data Histories 


RPM 

HERTZ 

GAS TEMP 

579. 

9.65 

- 

915.0 

15.25 

A2 

1246 . 8 

20.78 

A6 

1512.7 

25.21 

A10 

2434.8 

40.58 

A14 

3452.0 

57.53 

A18 


PRESSURE 

T 

X WALL 

S/ALL 

- 

A1 

A1 

A3 

A4 

A5 

A7 

A8 

A9 

All 

A12 

A13 

A15 

A16 

A17 

A19 

A20 

A21 


A comparison of most of the cycle -averaged surface tempera- 
ture histories and the corresponding surface heat transfer 


histories is presented in Fig. 13. 
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GAS TEMPERATURE 
EXPERIMENTAL OBSERVATIONS 

A sampling of the recorded gas temperature histories is 
presented in the histories of Appendix A (specifically 
Figs. A2, A6, A10, A14 , and A18). The noise level (both 
60 Hz and high frequency) was so low as to be negligible, 
except for two brief segments of the cycle when the syn- 
chronization breaker points actuated, at about -105 and 
+26 crankangle degrees, for four degree intervals. There 
was a moderate amount of cycle-to-cycle variation in the 
bulk gas temperature history, due probably to departures 
from one-dimensionality and random gas motion due to rock- 
ing of the piston and the induction of makeup air. The 
cyclic variation occurred primarily near top dead center 
(TDC) and 80 degrees to either side of bottom dead center 
(BDC) , and decreased somewhat with increasing frequency. 

The shape of the gas temperature history was essentially 
unchanged as driving frequency increased, except for an 
increase in amplitude and a slight shift in the crankangle 
at which the history maximized. Fig. 12 details this shift 
in bulk gas temperature maximum as recorded by the resis- 
tance thermometer: at 10 Kz the peak is at about -2.5 
crankangle degrees; at 60 Hz the peak is at -1.5 degrees. 






TEMPERATURE, scale is i.o f/unit 
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The bulk gas temperature swing increases over this fre- 
quency range from about 560 P° to about 630 F° . All gas 
temperature measurements were made with 3.8l u tungsten 
resistance wires. 

PEAKING OP GAS TEMPERATURE 

According to Pig. 12, at low driving frequencies the in- 
stantaneous gas temperature as measured by the resistance 
wire leads both the gas pressure maximum (-2.5 crankangle 
degrees) and the gas volume minimum, but at high driving 
frequencies the wire temperature peak lags the pressure 
peak slightly. 

The pressure transducer and its charge amplifier have 
very high frequency response: the Kistler model 566 elec- 

trostatic charge amplifier has a nominal 150 kHz upper 
limit; the 601H transducer has a flat response to within 
1 % up to at least 15 kHz (about the 300th harmonic of the 
highest frequency studied in this work) . The recorded 
pressure history can be considered as a faithful repro- 
duction of the chamber pressure. The resistance thermo- 
meter time constant is not accurately known, however, and 
some possibility exists that at high frequencies the tem- 
perature of the small wire may lag the gas temperature by 
a few crankangle degrees. Since both the gas pressure and 
the gas temperature are bulk averaged quantities, a simple 
first-law analysis may shed some light on how the wire tern- 


perature relates to the gas temperature. 
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The first law for a simple closed 
system is: 

PV = mRT . 

We can differentiate with re- 
spect to time and assume that 
the specific gas constant R is 
indeed constant : 



PdV V dP 

R dt R dt 


d_ 

dt 


(m-T) 


(19A) 


If we re- 
strict our- 
selves to 
the entire 
compression 
stroke, the 
dV/dt term 
is always 
negative. 

If we further restrict ourselves to the short time interval 
"B" between and including the pressure peak (about 2.5 de- 
grees before TDC) and the volume minimum (TDC, or 0 crank- 
angle degrees) in the compression stroke, the dP/at term 
is also negative or at most zero. Then in the domain "B", 
the quantity d(m‘T)/dt must be negative. For a constant 
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mass system, the gas temperature must be decreasing; 
then, the gas temperature must peak in the domain "A" 
prior to the pressure peak. 

If the system is complicated by the possibility of mass 
loss through blowby, the argument is not as simple. Exa- 
mine the instant at which the pressure reaches its maximum 

(i.e., -2.5°). The term V dP(-2.5°) is zero. The term 

R dt 

P dV(-2.5°) 

R cHF is » for fc yP ical conditions encountered in 

this investigation, approximately: 

P dV(-2.5°) y (388 psi )(-.0032 in^/cr. degree) 

R at ~ “ (53.3 ft lb f /lb m °R) 

P dV(-2 . 5° ) •v# 

r ^ “ -.00200 lb m °R/cr. degree. 

The equation now reads, 

-20 x 10 _ 1 , lb m °R/degree = T ^ + m ' (19B) 

The approximate size of the first term on the right side 
is, from the SIMUL simulation program, negative but small: 

T 2 * _d 290 °R ) ( .315 x 10 6 ^ ) 

cr. degree 

rn dm(-2 . 5° ) y -.0004 lb °R/cr. degree. 

1 dt " n 
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Then the term containing dT/dt in Eqn. 19B is negative 
when the pressure reaches its maximum for the case allowing 
blowby mass flow out of the system. 

The gas temperature, bulk averaged, is thus seen to peak 
before the pressure whether or not there is mass loss 
from the system. The experimental observation illustrated 
in Pig. 12 can be interpreted as follows: at high fre- 

quencies, the resistance wire temperature lags the gas 
temperature by two or three crankangle degrees . 


GAS PRESSURE 

Pigs. A3, A7, All, A15, and A19 present gas pressure his- 
tories at frequencies from 15 to 58 Hz. No noise was ever 
identified on a pressure trace. There was very little 
cycle-to-cycle variation; what occurred took place pri- 
marily in the lo w pressure period of the cycle- and perhaps 
was related to makeup air addition. As driving frequency 
increased, the cyclic variation in pressure history 
gradually decreased. 

The shape of the gas pressure history was the least affect- 
ed by increased driving frequency of all the histories 
measured. The mean pressure swing changed from 365.8 psi 
at 15 Hz to 36 k .2 psi at 58 Hz. The pressure maximum re- 
mained at -2.5 crankangle degrees with a maximum cyclic 
variation of perhaps one degree. 


SURFACE TEMPERATURE 
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Figs. Al, A4, A8, A12, Al6, and A20 present wall surface 
temperature histories measured at a point near the center 
of symmetry of the head at frequencies from 9 to 58 Hz. 

A summary of most of the surface temperature histories 
(each history being the average of several cycles) with 
the corresponding heat transfer histories is given in 
Fig. 13. At low speeds (up to about 15 Hz) the wall sur- 
face temperature peaks unexpectedly late. The lag of the 
surface temperature peak (with respect to a fixed reference 
such as TDC, or the pressure history peak) increases with 
increasing speed. At high speeds (35 Hz to the maximum of 
the equipment, 58 Hz), the wall surface temperature peak 
is at or very near to top dead center. At intermediate 
frequencies the two trends overlap and a clear double 
peaking is observed in head (or wall) surface temperature. 

A plot showing the trends of the crankangle position of 
this temperature peak with frequency is shown in Fig. 14. 

An interpretation of these data is deferred until later. 
These histories represent very small signal levels: a 2 F° 

swing in surface temperature produces a voltage swing of 
only 60 microvolts. As a consequence, noise suppression 
was important; the steps taken to increase the signal-to- 
noise ratio were discussed earlier. The observed 60 Hz 
ripple was less than 5% of the signal size, and the high 
frequency noise level was held to less than about Q% of 
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There was considerable cycle-to-cycle variation in the 
surface temperature data; for this reason, several cycles 
were recorded at each speed and the results later were 
averaged for interpretation. In spite of cyclic variation 
in the surface temperature histories, very definite 
patterns emerged from the data; these patterns are dis-r 
cussed later. 

SURFACE HEAT TRANSFER 

LISTING OF DATA 

Surface heat transfer histories are given in Figs. Al, A5, 
A9, A13, A17, and A21. The method of calculation of these 
histories from the surface temperature histories is ex- 
plained in the next section. In spite of the -cycle-to- 
cycle variation in the surface temperature data, the sur- 
face heat transfer histories averaged over a number of 
cycles show definite patterns of change with increasing 
driving frequency. This is best seen in the summary figure. 
Fig* 13, together with the very low speed data of Fig. Al. 

COMPUTATION OF INSTANTANEOUS HEAT TRANSFER 


Statement of Method 


The instantaneous heat transfer at the compression 
chamber surface was obtained by the application of one- 
dimensional semi-infinite solid theory to the mea- 
sured chamber surface temperature. It would, in principle 
be possible to measure instantaneous surface heat transfer 
rates at the temperature levels and cyclic frequencies 
covered in this work with recently developed heat flux 
transducers*, but this course of action was not taken for 
a number of reasons: 


1) xt is possible to calculate the required heat 
transfers accurately, though at the expense of 
time and effort, from the measured surface tem- 
perature history. 


2) Available transducers are on the order of 1/4 
to 1/2 inch in diameter, and instrument spacing 
would mean that the temperature and the heat 
flux were measured perhaps one inch apart in a 
2.3125 inch chamber. 


3) Available transducers usually have ceramic or 
pyrex backup materials, severely disturbing the 
one— dimensional heat flow and making dubious a 
correlation between heat flux and surface tem- 
perature . 


4) Instantaneous heat flux recording would require 
another channel in the already filled display 
(a four-channel oscilloscope). 


For these reasons and others, the calculation route to 
instantaneous surface heat transfer was taken. 


*Heat Technology Laboratories in List of Manufacturers 
and Addresses 
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The basic principle of calculation is as follows: the 

surface temperature history is taken from a film record 
and converted to an array of ordinates and associated 
crankangle times. The history is taken as being given 
by these points and linear line segments between the points. 
A Fourier analysis is performed to break the surface tem- 
perature history down into its harmonic elements; a limit 
for the number of harmonics is chosen (a number less than 
the number of original ordinates); the higher harmonics 

are not used. Semi-infinite solid theory is applied to this 
harmonic synthesis, and an expression for the instantaneous 
temperature profile into the wall is generated. Instan- 
taneous surface heat transfer is calculated from a differ- 
entiation of the temperature profile, evaluated at the 
plate (i.e., head) surface. 

This technique is the same as used by Overbye (33) and 
others. It involves fairly standard heat transfer philo- 
sophy; there are, however, significant differences in cal- 
culation and application between the method used and that 
of Overbye, resulting in greater flexibility and accuracy. 
Credit for the origin of these differences and their trans- 
formation into a method of calculation should be given to 
T. LeFeuvre*. A more detailed study of this method will be 
found in LeFeuvre's Ph.D. dissertation (1968). 

*Ph.D. candidate. Department of Mechanical Engineering, 
University of Wisconsin 


One-Dimensionality Approximation 
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Heat transfer is assumed to be one-dimensional, normal to 
the plane of the head surface, throughout this work. The 
surface thermocouple from whose signal the flux is calcu- 
lated is close to the center of a half-inch thick plate of 
cast iron which constitutes the head. Fig. 15 shows that 
the resistance thermometer probe is the only ijearby object 
which could disturb the heat flow (other than the instru- 
ment plug-head thread interface), and it is four thermo- 
couple diameters away. The washer-shaped thermocouple 
junction is backed by iron, the same material as the head. 
It is essential* that the thermocouple backing material 
be identical to that of the head if the output signal is 
to be representative of the typical surface temperature of 
the head. The central oxidized nickel or constantan wire 
does constitute a different thermal resistance in parallel 
with the iron resistance ( k constantan = 17 Btu/hrft°F, 
k iron = 29 Bt u/hrft°F), but the wire area normal to the 
flux direction is only 1.355 of the junction area. 

Semi-Infinite Solid Approximation 

The use of semi-infinite solid laws is valid if the. surface 
temperature fluctuations are damped to negligible ampli- 
tude before reaching the exterior surface. Overbye de- 


*Periodical #TB— l6l, Nanmac Corporation 


Two Views of Underside of Head 

Instrumented with resistance wire at top and one 
surface thermocouple on the right in central plug, 
and a Kistler pressure pickup in the plug at right 
in one of the semicircular volumes in the head. 




Figure 15 - Underside of Head 
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tected no perturbations in exterior temperature at low 
frequency operation (low frequency harmonics penetrate 
deeper than higher harmonics) with a 35 P° swing on the 
inside of a cast iron wall less than 3/8 inch thick. This 
experimenter made the same observation with better elec- 
tronics, a 2 F° inside swing, and a half-inch wall. 

According to Schneider (67), the envelope of a temperature 
history which is periodic at the surface decays with depth 
into a semi-infinite solid according to: 

Envelope size ^ exp (-x/w/2a ) 

where w is the frequency of the lowest harmonic, 

a is the thermal diffusivity of the material, and 
x is the distance into the material. 

If we set a 1000:1 diminution as the standard for the semi- 
infinite solid, the critical depth for "complete” damping 
Is given by x in the following expression: 

V 

exp (-0 Vw/2a) ■ 1000*exp ( -x^/w/201) , or. 


x = +/2a/o log 1000. 

C 6 


( 20 ) 


For cast iron (a ■ 0.666 ft 2 /hr), several values of x 

c 

for various gj are summarized below: 


w £r£m)_: 
x c (in.): 


_800 

.174 


1000 

.156 


1500 

.127 


2000 

.110 


i°oo 

.078 
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It Is clear that semi-infinite solid theory can be applied 
to a half-inch slab of cast iron with confidence. 

Theory 


Fourier Series Approximation 

A wall which can be taken as a semi-infinite solid sees 
a steady periodic temperature history at its surface. We 
assume that the oscillation has been maintained long enough 
so that the original transient state due to starting the 
cyclic surface temperature has built up into a steady type 
of periodic flow. 


The first concern is to determine the time-dependent tem- 
perature profile in the wall. The heat transfer litera- 
ture contains many derivations for temperature profiles 
induced in slabs and cylinders by surface temperature var- 
iations expressed in terms of elementary functions. In 
the case of experimental periodic temperature variations 
not easily expressed in terms of simple functions, the best 
procedure is to reduce the experimental curve to an equiva- 
lent trigonometric series. The profile into the wall for 
each series term is obtained by application of the Fourier 
heat equation: 


3_t 
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a 


aft 

3x 2 


( 21 ) 
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where: t is temperature, °P, 

0 is time, hours, _ 

a is thermal diffusivity, ftVhr., and 

x is distance into the wall from the surface, ft. 


Let us consider an experimental temperature-time history at 

the metal surface consisting of N+l temperature ordinates 

t (0,9 ) at corresponding times 0 . The abscissas 9 need 

n n 

not be equally spaced. 


We take the function t(0,e) to be defined by the array of 
N+l points and linear line segments between the points; 
for reasons not discussed here, the ordinate of the first 
point is discarded and the value for the N+lst point is 
substituted so that the end point ordinates are equal. A 
Fourier analysis will provide a trigonometric equivalent 
to the function; we will limit ourselves to the first K 
terms : 


K 

t(O,0) = A qw + ' [A k cos(ku>9) + sin(kw0)] (22) 

k=l 


where t(O,0) is the approximated temperature function, 
k is the harmonic number, k = 1,2,3, ...K, 
w is the frequency in radians/hour, 

0 is the time in hours, 

A qw is the average of the given ordinates, 

A k is the amplitude of the kth cosine harmonic, 

B k is the amplitude of the kth sine harmonic. 


The A and B coefficients are given by: 
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ow 


1 

0 


N 


t(O,0) d0 


u 0. 


oe 


N 


2 

0 


t(O,0) cos(0irk/K) d0 , and 


0 


(23) 


(24) 


B k ** 0 t(O>0) sin(0*k/K) d0 (25) 

^0 

where k = 1,2,3, 4, ...K 
0 O ^ 6 * 

0 E ^ 0 N “ 0 o^i ( 0 q would generally be taken = 0). 


A more convenient, but less accurate, definition for the 
coefficients is more often used in which the integral 
definitions for A QW , A^, and become sums over the N 
ordinates : 

N 

A «w = ir 2Z V^V • (26) 

n=l 

N 

A k 5 I t n (O » 0 n ) cos (n2Trk/K), (27) 

n=l 
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N 

t n (0 »°n ) Sin (n 2 ^/K). (28) 

n=l 

In this form, however, the function is broken into a series 
of steps; the height of each step is the ordinate value at 
the end of the step, and the trigonometric functions are 
evaluated at the step ends. The approach used. here re- 
places the steps with linear ramps and actually performs 
the integrals. 

Applying the separation of variables technique suggested 
by Schneider (6 7), with the boundary conditions: 

X 

t(0,6) = A qw + ^ ' A^ cos(kwQ) + sin(ku>0), (29) 

k=l 

and t (L, 0) = A . 

aj 



where L is the wall thickness in feet, and 

A aj . is the external (air jacket) temperature, °F, 

we get the steady state solution. 


t ss (x,9) = A ow - l Aa h ., ! £ 


ow 


(30) 


and the transient solution, 

K 

\ I \ / 

t tr (x,6) = \ exp(-px) [A k cos(ko)0-qx) + B k sin(kw0-qx) ] . 

k=l 


(3D 


when the p and q are evaluated by subjecting thi3 
solution to the test of the Fourier heat equation, we 
find that 

p = q = /ka)/2a (32) 

and the complete solution for the temperature profile 
induced in the wall is: 


t(x,6) 


ow 



+ 


( 33 ) 


+ 





cos(kw0 - x 



[ sin(kw0 




This equation satisfies the boundary conditions at x *» 0 
and x ® L: when x = 0 it reduces to 


K 

t(O,0) - A qw + ' [A k cos(kwO) + B k sin(k(D 0 )] ( 34 ) 

k=l 


and when x * L, recalling the infinite plate assumption 
that the time dependent portion of the solution is negli- 
gibly small, 

t (L, 0) = A aj . (35) 

Differentiating the equation for the induced temperature 
profile with respect to depth into the wall, we get: 
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3t(x, 0) 

dX 


^ow ~ ^a.i 


K 



ka) 

2a 


exp 


~Jg 


k=l 


(A k - B k ) sin(kw0 - x/§£) + 


(>B k -A k ) cos (kw0 - x/i) 


k k 

Setting x equal to zero. 




(36) 


3t(x ,'Q) 


3x 


A ow ~ A a,j 


x=0 



(Bk-Ak)sin(ka)0) + 
(Bk+Ak)c°s(kw0 ) 


(37) 


Since the surface heat transfer is defined from the metal 
side of the surface by 


q 


surface 


= - (cond.). 8t( ;' 9) 


3x 


x=0 


(38) 


the heat transfer rate, in units of Btu/hr ft , is given by 


^surface " cond *l 


^ow ~ ^a,1 



K 


+ /Sr • >Jk 



k=l 


(B k -A k )sin(kw0) + 
(B k +A k )cos(lcw0) 


(39) 


The first term of this expression represents the steady- 
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state heat transfer rate and the series summation 
constitutes the transient heat transfer rate at the wall 
surface . 

Converged Fourier Series Approach 

When approximating a periodic function with a Fourier series 
with a "reasonable" number of terms, a series .representation 
often "overshoots" any steep-fronted segments and oscillates 
about the function. Hamming (68) points out, for example, 
that when approximating a square wave of unity amplitude, a 
ten-harmonic Fourier fit overshoots by nearly 9 %, then 



oscillates about as shown in the sketch above. This be- 
havior is often termed the "Gibbs Phenomenon . " 
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In treating experimental wall temperatures to Fourier 
analysis with 50 harmonics, and computing instantaneous 
surface heat transfer, the resultant histories appeared 
to contain oscillations of regular frequency as shown in 
the lower half of Fig. 16 . These oscillations were attri- 
buted to the Gibbs Phenomenon, as increasing the number of 
terms in the analysis from 50 to 70 terms decreased their 
amplitude. 

A number of mathematical "tricks" exist to smooth the 
Fourier representation, generally at the expense of high 
frequency response. If each coefficient A. and B. is mul- 

K 

tiplied by the so-called Lanczos Smoothing Factor, 


°k 


sin 7 rk/K 
uk/K 


( 40 ) 


where k » 1 , 2 , 3 , ...K, 

the first overshoot in the previous sketch is reduced from 
0.08949 to 0.01187, and the first minimum is reduced from 

0. 04859 to 0.00473 (for the 10 -harmonic representation, 

1. e., K = 10). The factor is very close to unity for the 
first harmonic and decreases to zero for the Kth harmonic, 
as' seen in Fig. 17 . 


Applying the Lanczos factors to the harmonic coefficients 
in the calculation of the wall surface heat flux damped 
the regular oscillations, as shown in the top half of Fig. 
16. Essentially the same result could have been achieved 





NIC NUM3E 
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by increasing the number of terms K to over 100, but at 
the expense of increased computer time. For want of a 
better name, this computation technique is herein termed 
the "converged Fourier" method as contrasted to the 
"Fourier" method. 


With the application of this factor, the surface tempera- 
ture history becomes, instead of Eqn. 29, 


*«••> ■ A OW + 


K 

' ~Sk/K^ j^ A k cos ( ka)Q ) + B k sin(kw6)J- 

k*l 

(41) 


and the temperature profile into the wall is now: 


t(x,e) = a ow 


( ^ow ~ A a,j ) x 
L 


(42) 



sin(7rk/K) 

irk/K 


A, cos(kw8-XiV ^) + 

K v 2a 

B k sin(kwe-xv^) 


iiie instantaneous heat transfer rate at the surface can be 
evaluated from: 


‘converged, surface 


cond. 


-A. 


ow 


- A 


Si 


(43) 


52 




A^) sin(kwe) 


(B k + A^) cos(koiO) 



As the number of terms considered, K, increases, the small 
differences between the measured wall surface temperature, 
the "Fourier" fit (Eqn. 29), and the Lanczos smoothed or 
"converged" Fourier fit (Eqn. 41) will decreas-e. For 50 
terms, the average error of the fits with respect to the 
measured curve is on the order of 1 % of the total swing 
(i.e., about .02 F°). 


CALCULATED FILM COEFFICIENT 

Newton's Law of Cooling, expressing the rate of heat trans- 
fer across a film at the interface of a fluid and a solid. 


q/A 


surface 


h-(T„ 


luid 


T 

solid 


surface 


)*, 


( 44 ) 


has been commonly applied for many years in the study of 
free and forced convection. Its use is quite conventional 
and proper in many practical steady-state applications. 
However, caution must be observed in using this equation 
to describe the heat transfer rate in a system with rapidly 
changing conditions. The heat flux q/A is evaluated at the 
solid surface. When the properties of the film next to the 
surface are such that a phase lag occurs in the response of 
the wall surface temperature with respect to a gas tempera- 
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ture change, the film coefficient h can assume negative, 
zero, and infinite values. This effect is described 
briefly by Overbye(33) .. A more recent and complete ex- 
position on negative and infinite values of h is given in 
a text by Y. V. Kudryavtsev ( 69 ). The theoretical models 
in this text incorporating capacitive energy storage in 
the boundary layer exhibit this behavior (see Pig. 58 for 
the RC Model, explained in Section VIII. B.. 4. ) . 

The experimental data, too, produce the peculiarly be- 
having h. Pig. 18 depicts the behavior of the coeffi- 
cient calculated from: 


experimental (T 


_a/A_ 


- T T 

gas wall surface ' 


(45) 


When h is positive, the heat flux and temperature differ- 
ence agree in sign. When h is negative, they differ in 
sign. When h is zero, the heat flux is zero. Wh en h is 
infinite, the film temperature difference is zero. 


The coefficient h assumes infinite values whenever the 
surface heat flux does not change sign precisely when the 
temperature difference across the film changes sign. 

AERODYNAMIC COMPLICATIONS 

Ideally, all gas motion and property gradients in this 

s. 

investigation take place in only one dimension, normal to 
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the head surface. The theoretical studies are all based 
on models in which a boundary layer is bounded on one side 
by a wall (the head) and on the other side by a spatially 
homogeneous gas mass (the bulk gas) in which there are no 
temperature gradients. The resistance thermometer, pro- 
truding about 0.100 inches (about 3 mm) into the gas, mea- 
sures this bulk gas temperature. When several resistance 
thermometers were inserted the same depth into the gas at 
different places, the temperature histories seen were iden- 
tical, showing the independence of temperature from r and 6, 
However, in order to have a gas mass with uniform density, 
there must be some small-scale gas motion in the bulk to 
mix the cool gases near the cylinder walls with the warmer 
"core" gases. It was therefore thought desirable to mea- 
sure i/he head surface temperature at a number of positions 
simultaneously to see if the temperature history varied 
appreciably with location on the head surface.' To this 
end, a head was constructed with nine surface thermo- 
couples distributed as shown in Fig. 19. Surface tem- 
perature data were taken over the frequency range for 
three conditions as described below. 

CYLINDRICAL CHAMBER 

The head of Fig. 19 differs from the original head (see 
Fig. 2) in that the two semi-circular clearance volumes, 
or "ears" (one to accomodate the pressure transducer, one 


Underside of modified head with nine 
surface thermocouples, after running 
for 2 hours at 2000 rpm. The black 
specks are teflon residue from rings 
and cylinder walls. 


Fig. 19 
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on the opposite side of the chamber for symmetry), were 
not needed, and therefore not included. As a consequence, 
this compression chamber was truly cylindrical and the 
fluid did not experience any "squish" flow (into the ears 
during compression, out of the ears during expansion). The 
nine surface thermocouples all experienced essentially the 
same histories as described by Figs. Al, A 4 , a 8, A12, Al6, 
and A20 for driving frequencies up to 22 Hz. Above this 
frequency the history was much the same as in the figures 
except that a "parasitic" oscillation at about the fifth 
or sixth harmonic of the fundamental appeared superimposed 
upon the history, growing in amplitude with speed until at 
58 Hz its swing was about 4 0 % of the main temperature 
swing. Fig. 20 shows the crankangle locations of the 

peaks of the "primary" and the "parasitic" temperatures on | 

i 

the head surface; this figure can be compared directly to 
Fig. 14 for the original head. Fig. 21 presents an example 
of the surface temperature behavior at 3175 rpm (53 Hz), 
showing the superimposed oscillation. The "parasitic" 
wall surface temperature oscillation is presumed to be due 
to changes in the heat transfer coefficient induced by 
cyclic air flow pattern in the chamber. All nine thermo- 
couples reported this temperature oscillation, although 
its amplitude varied somewhat about the head surface. 

I 

i 

’ i 

r 

\ 

i 

! 
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CRANaANGLE 


CYLINDRICAL CHAMBER WITH CLOSED MAKEUP PORTS 
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On the suspicion that air flow through the small makeup 
fluid ports in the cylinder wall bottom was serving as a 
driving mechanism for the recurring air flow (which was in 
turn inducing the periodic "parasitic" heat transfer and 
surface temperature swing), the holes were plugged and a 
set of data runs made. The size of the temperature and 
pressure swings decreased with time as blowby past the 
piston rings caused the system mass to decrease, but the 
wall surface temperature oscillations were not decreased 
in ratio to the overall wall surface temperature swing. 

The cyclic air flow is therefore not totally dependent on 
makeup air flow for its driving function. An alternative 
driving mechanism is the scraping of the cool boundary 
layer off the cylinder walls by the piston and the mixing 
of this denser mass with the bulk gas. 

CYLINDRICAL CHAMBER WITH SCREENING 

A screen of fine-mesh nylon was constructed 0.120 inches 
from the head interior surface and parallel to it to break 
up the cyclic air flow patterns and reduce their effect on 
surface temperature and surface heat transfer. This con- 
struction is illustrated in Pig. 22. The screen had two 
layers: a support of Nitex 35-500 nylon monofilament (32 

filaments/inch, 300 micron filament diameter, 500 micron 


SURFACE THERMOCOUPLES 



ASSEMBLY WITH SCREEN IN POSITION 
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opening, 392 open area), and a fine mesh of soft nylon 
screening of about 100 micron opening, 150 filaments/inch. 
The screen was fastened, around the circumference of the 
chamber and was supported at five locations on the face of 
the head by small plastic tubes (see Pig. 22). 

The parasitic oscillations on the surface temperature his- 
tories were reduced by 50 to Qu% in amplitude ‘from the 
previous runs with the unshielded "cylindrical" chamber. 
The maxima of the small oscillations occurred at the same 
crankangle locations as before, as indicated by Pig. 23. 
This figure can be compared directly with Pig. 20 for the ' 
unshielded chamber and Pig. lH for the original head. 

CONCLUSIONS ON AIR MOTION 

The three data sets discussed above shared with the ori- 
ginal set of Pigs. Al, A*», A8, A12, Alo, and A20 (or, in 
summary. Pig. 12) the same general trends in head surface 
temperature history: a lagging surface temperature hump 

at low driving frequency; a peak near piston top dead cen- 
ter (0 crankangle degrees) at all high speeds; a double 
peak combining the two phenomena at intermediate speeds. 
The original set was completely free of the "aerodynamic" 
effects of the last three, and was taken to represent 
simple one-dimensional flow normal to the head surface. 

The oscillations were probably due to cyclic air flow 
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patterns sustained by very small pressure gradients 
in the bulk gas; the "squish" effect of the small semi- 
circular volumes (0.98$. of the total clearance volume) 
was sufficient to break up these resonant air flow patterns 
in the original tests. 


X 
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VIII. THEORETICAL MODELS 

A number of theoretical models were developed to treat 
several domains in the system. All of the models take the 
form of digital computer programs, and each employs some 
type of relaxation or cyclic iteration technique to proceed 
from a guessed initial state to a cyclic history of states 
which satisfies the pertinent energy and continuity equa- 
tions . 

One of the theoretical models, termed the Adiabatic Plane 
Model, was considerably superior to the other models; it 
will be discussed in detail on the following pages. The 
five secondary models, which preceded the other chronologi- 
cally, will be briefly described at this point. A more 
complete listing of the assumptions, methods of calculation, 
and results of the secondary models can be found in the 
Appendix section of the author’s Ph.D. thesis.* 

The first of the "secondary” models is named SIMUL, and 
treats the entire working fluid and the surrounding hard- 
ware: head, piston, and cylinder walls. This model was 

successful in predicting the steady periodic behavior of 
the working fluid treated as a homogeneous mass: bulk gas 

temperature and pressure histories, and heat flux histories 
through the surrounding metal surfaces (based on constant 

*Ph.D, thesis, Daniel W. Wendland, U. of Wis., Jan. 1968 
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internal metal surface temperatures) were calculated. 

The remaining four "secondary" models consider the head 
(sometimes termed the "wall") and the boundary layer ad- 
jacent to it on the compression chamber, or working fluid, 
side. In each, the intent is to examine the combination 
of various energy terms in the boundary layer with the 
energy flux entering the external edge of the boundary layer. 
The metal surface heat transfer rate is the result of this 
interaction of energy fluxes, and the metal surface temper- 
ature responds to changes in this heat flux. Each of these 
models was successful, to a varying degree, in predicting 
the metal surface temperature and surface heat transfer 
responses to gas temperature and pressure driving functions. 
They are not presented here, however, since the Adiabatic 
Plane Model rests on a far surer theoretical base and makes 
a better prediction of the metal surface thermal response. 
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ADIABATIC PLANE MODEL 

In each of the secondary models of either the entire work- 
ing fluid and the surrounding metal walls or of the head 
and its gas-side boundary layer, a considerable amount of 
empirical information was necessary to make the analysis 
possible. This empirical information was typically in the 
form of a time-varying spatial or energy flux boundary 
condition. In the total system model (SIMUL) , the heat 
transfer rates at the metal surfaces were taken to be pro- 
portional to the Eichelberg heat transfer coefficient and 
to the bulk gas - wall surface temperature difference, 
both time-varying quantities. In the four models of the 
boundary layer - wall system, the Eichelberg coefficient 
was again used, and the time-varying boundary layer thick- 
ness was estimated around the cycle. The latter was ac- 
complished by specifying the value of the Nusselt Number 
based on the instantaneous heat transfer coefficient and 
the instantaneous boundary layer thickness; this value, 

NUg, was generally held constant around the system cycle. 

It was occasionally necessary to assume that the velocity 
profiles from piston to head were linear. Since some of the 
assumptions do not have strong theoretical bases, it was 
decided to construct a model for the system which draws 
its resources from basic thermodynamics and contains a 
bare minimum of simplifying assumptions and "engineering" 
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data. This theoretical construct is known as the Adia- 
batic Plane Model. 

The system is the working fluid constrained by the head, 
cylinder walls, and the moving piston. The first major 
assumption, 

1) the system is one dimensional, normal to head 
and piston, 

ignores heat transfer to the cylinder walls and its in- 
fluence on the gas properties, and can be defended on two 
grounds: first, the experimental system to which oh is 

model is an analogue was designed so that the gas motion 
would be as close to one-dimensional as possible; second, 
our interest is not so much in a gas being compressed and 
expanded in a particular piston-cylinder combination, but 
rather in the unsteady thermodynamic state of a gas ex- 
periencing pressure and temperature fluctuations and ex- 
changing energy with a plane solid wall. 

The second assumption is that: 

2) the head surface temperature is constant. 

The constancy of the metal surface temperature is not too 
bad an assumption, since the expected surface temperature 
swing is less than 1 % of the gas temperature swing. The 
assumption of constancy is not fundamental in the analysis, 
in any case, and it will be shown later how a known varia— 


ble surface temperature could easily be introduced to 
the system model. 

The third major assumption is that: 

3) ' there is a plane midway between piston and head 
at which the gradient of gas temperature with 
distance is zero. 

This would follow from symmetry if we made the piston 
and head surface temperatures constant ana equal. How- 
ever, even if the piston surface temperature is not quite 
equal to the head temperature (the simulation program pre- 
dicts it to be some 50 to 80 P° higher), its time variation 
is on the order of a few degrees — less than 1 % of the gas 
temperature swing. The time-varying location of the adia- 
batic plane (or plane of zero temperature gradient) may be 
slightly different from the instantaneous position of the 
plane constrained to be halfway from piston to head. How- 
ever, the temperature gradients in this region are very 
small (see figures in this section) and the error intro- 
duced by this boundary condition is small. 

It is simplest, then, (and least expensive in terms of com- 
puter time) to consider the system to be bounded by the 
head and by a plane moving midway between head and piston. 
At the plane, the heat flux is zero; at the wall, the tem- 
perature is constant . 


The final major assumption is that: 
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*0 the fluid pressure is a function solely 
of time. 

■<* 

This assumption is quite reasonable at cyclic frequencies 
^ fsw hundred Hz, since the rate of change of chamber 
pressure is relatively slow, and the system can quickly 
•'smooth out" any small pressure gradients existing in the 
gas. As frequency increases, fluid velocities begin to 
approach sonic values, and pressure gradients can stand in 
the gas . 

c 

Other assumptions are as follows : 


5) 

6 ) 


J -he fluid is a periect gas with the molecular 
weight and basic properties of air. 


lhe gas specific heat and thermal conductivities 
can be approximated by linear temperature func- 
tions: at temperatures below 500°F, 


k gas 

c 

p 

*gas 
and at 


.0132 + (2.36 x 10“ 5 ).T (Btu/hr ft °?) 
: .2400 + (1.71 x 10“ 5 )*(T-90) (Btu/lb 

temperatures above 500°F, 



k 

gas 

c 

p gas 


0250 + (1.70 x 10“ 5 ) (T-500) 
.2470 + (3.20 x 10" 5 ) (T-500), 


where T is the local instantaneous gas temperature 
in degrees P. 


7) The mass of the system is constant. 

8) The fluid is inviscid. 

9) ±iiQ fluid is not subject to inertia forces; the 
kinetic energy is small. 


10) There is no gravity field, or, the one dimension 

of the system is perpendicular to the direction in 
which gravity is acting. 


71 


Assumption 7 becomes important when one compares the 
results of this model with the experimental data, taken in 
a working fluid whose mass changed slightly around the 
cycle. The ninth and tenth assumptions state that the 
energy dissipation due to fluid friction and the fluid's 
kinetic energy are negligibly small. 

THE MATHEMATICAL FORMULATION 

The gas, with total mass M, is subdivided into N elements, 
or laminae, as shown in Fig. 2H. Each lamina has a con- 
stant mass, one Nth of the total mass, and is at a time— . 
dependent temperature T throughout. All laminae share 
the same pressure history by assumption 3. 

There are N+l equations describing the thermodynamic state 
of tnis system at any time t. The first equation says 
simply that the sum of the volumes of the laminae at the 
time t must be equal to the volume contained at that time 
between head and plane (i.e., half the piston displacement 
volume at time t). Since the system is one-dimensional, 
all equations are written per unit area. The area is taken 
as unity for convenience. Then the first governing equa- 
tion is: 

N 

Tv*) - V*> (46) 


n=l 



7iC. 2b 


The Adi 
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where V (t) is the volume of the nth lamina at time 
t, and 

V, (t) is the volume between head and plane (in 

u -3 p 

units of inches ^/inches ) at time t. 

The other N governing equations are energy rate balances 

at time t for the N laminae. For the nth lamina, with 

thickness Ax^ (equal to V , since the units are inches^/ 

2 

inches ), the simplified first law for the general element 
n at time t is : 


l n 


in 


00 - q„ 00 ■ 


‘n 


me 


out 


dT (t) 
1 1 

dt 


+ P(t) 


dV (1 
n 

dt 


47) 


where m is the mass of element n per unit area. The units 
of this equation are energy/length^ time. 


The last term, which expresses the rate of work, can be 
changed with the help of the perfect gas law: 


PCt) 


av n (t) 

dt 


PCt) 




pityj 


(48) 


P(t) 


dt 


- mRT n (t) 


T 
L n 


1 

TtT 


dT n (t) 


1 dP(t)1 

P(t) dt J (49) 


Applying the relationship for a perfect gas, c p = c y + R, 
the first law for the nth element is: 


v <*> - % 

in . out 


dT (t) 

me • 


mRT ( t ) . 

• n v ' cP(t) 

4*1 Vi^ i 

p dt 

md 

! 


_"pTtl dt 


( 50 ) 
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The first term in this equation represents the energy 
input rate by conduction to the nth mas 3 , or lamina, as it 
and the observer moves in the system; the second is the 
output rate. The third term is the rate of change of en- 
thalpy in the nth chunk, and the last is a work-related 
term. 

THE FINITE DIFFERENCE EXPRESSIONS 
Development of Equations 

In the finite difference versions of the constraining 
equations, the choice of forward difference, backward 
difference, and central difference expressions for the 
first derivatives in the individual terms is guided by 
stability and convergence considerations. The forms se- 
lected below are the result of a series of trials with 
different forms for the terms which resulted in varying 
degrees of instability or stability. 

The work-related term is expressed as a central differ- 
ence. If we use j to represent the present time (at which 
all quantities are known or assigned), j-1 to be the 
immediate (and also known) past, and J+l to indicate the 
immediate (and unknown) future, then the term becomes: 

Work-Related Term 5 [ p(J+1)-P(J-1) 

P(j) L 2At 


.( 51 ) 


where At is the duration of the time interval. 
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The enthalpy term is expressed in forward difference form 
so that the temperature at time j+1 can be evaluated: 

Enthalpy Term = me 

P 

The use of the forward difference form is fairly standard 
practice in this type of evaluation. The central differ- 
ence form for this term was checked out, but its use re- 
sulted in rapid instability at all speeds considered. If 
the time interval is chosen to be sufficiently small, the 
error associated with the forward difference form for this 
term is small . 


T n (J+l) - T n (j) 
At 


( 52 ) 


The conduction rates are expressed for time j 


Q n U) = k 
n in 


T n+ i<J) - T n (j) 


Ax + ( j ) 


( 53 ) 


and 


q n (j) a k 
n out 


VJ> - Vi<J> 


Ax~( j) 


( 54 ) 


where Ax (j) is the distance from the center of element 
n+1 to the center of element n at time j, and Ax“(j) is 
the distance between centers of elements n and n-1. 

Because Ax n (j) and V n (j) are numerically equal (since depth 


is really volume per unit area), Ax + (j) can be written: 
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Ax + (j) 


Ax n+l ( J > , Ax n ( ^ 
2 T 2 


or, by the perfect gas law, 

Ax+( J> ■ 5K37 CT n + i<J> - 

and similarly. 


Ax"( j ) 


mR 

2P(J) 


CT n (j) + T n _ 1 ( J)]. 


Then the input conductive flux is: 



2_kP(J> T n+l ( ' I) " T n <J> 
mR L T n+l lJ > + V J > ’ 


and the output conductive flux is: 



(J) 


2kP(,j) 

mR 


- T n-l<J> 

VJ> + T n-l ( ^J ’ 


Combining into the complete finite difference form 
the nth, or general, element, 


2kP(,n WJ) - T n ( J> 
^ L T n+l U) + T n U) 


T n^ ~ Vi 
V^ + Vl 


(J) 

Til 


(55) v 


(56) 


(57) 


(58) 


(59) 

for 
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me 


T n (j+1) - T n (J) 

mKT n (j) 

P(j+1) - P(.1-l) 

At 

“ 

TTjD 

2At 


( 60 ) 


The analogous equation for the Nth, or last element re- 
flects the adiabatic wall boundary condition, q (J) * 0: 

n in 


2kP(.j) 

mR 


V-l) + 


V-> +1 > - t n (j) 

mRT N (j) 

P(J+1) - P(.1-l) 

At I 

P(J) 

2 At 

- _ 


.(61) 


The analogous equation for the first element reflects the 
constant wall temperature boundary condition in the follow 
ing way: since 


out 


U) 


T l ( J ) - T wall 


Ax“( j ) 


(62) 


and the appropriate Ax is the distance from the center of 
element 1 to the wall surface, or 


Ax"(j ) 


V 1 (j) mRT 1 (j) 

2 8 “TTCJT 


( 63 ) 


then the first law equation for the first element is: 


2kP(,j) 

mR 


T 2 (j) - T x (j) 

T 2 (j) + T x (j) " 


T X (J) - 


wall 
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There is nothing in these equations which says that T ,, 

wall 

must -stay constant; it was chosen to be so because the 
wall surface temperature history is unknown and the con- 
sequences of making it a constant are probably minor. If ’ 
a wall surface temperature history is available, it can 
easily be expressed as a subscripted variable and intro- 
duced to Eqn. 64. 

The N+l equations determining the thermodynamic behavior 
of the system are as follows: Eqn. 46, the volume con- 

straint; Eqn. 6l, the first law for element N; Eqn. 64, 
the first law for element 1, and the N-2 versions of Eqn. 

60 for the N-2 general elements. The unknowns are also 
N+l in number: the pressure P(j+1) at the immediate 

future time, and the future temperatures of the N elements, 
T n (j+1), n=l,2,3,...N. 

Two questions were encountered during the development of 
this model regarding: 1) the basic form of Eqn. 47 and 
2) whether the system is completely determined by the N+l 
equations listed above. These questions are assuaged in 
the following sections . 
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Displacement Work 

As the adiabatic plane mirrors the piston motion (but at 
"half scale"), the gas laminae are compressed and dis- 
placed. Their centers of gravity are moved, leading to 
the speculation that the work term in Eqn. ^7 should some- 
how contain an additional displacement or "pumping" work 
factor. 


The correct form for the work term consistent with the 
assumptions contained in this model can be deduced from 
a simple thought experiment. Consider a medium in one- 
dimensional space, in which the fluid pressure is a func- 
tion only of time. Planes 1 and 2 move with instantaneous 
velocities u^(t) and 


u 2 (t) in this fluid. 
The instantaneous 
rate of work on plane 
1, with area A, is: 


(System 1,2) 


u 1 (t) 


u 2 ( t ) 


( 1 ) 


P(t) 


( 2 ) 


W 1 (t) » P(t) • A*U 1 (t) . 


( 65 ) 


Similarly, the rate of work on plane 2 is: 


w 2 (t) • P(t ) -A*u 2 (t) . 


( 66 ) 


The net rate of work on the mass contained between the 


planes is the difference: 
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w 1>2 (t) ■ P(t) • A* [u 2 (t) - u 1 (t)] (67) 

Now, [u 2 (t) - u 1 (t)] is the relative velocity of the sides 
of the subsystem 1,2. Then the rate of change of volume 
of the subsystem could be expressed: 

dV (t) 

^ * A • [u 2 (t ) - u^t)] (68) 

and the net work rate for the lamina 1,2 is 

dV (t) 

w lj2 (t) = P(t) - , (69) 

which is precisely the form used in Eqn. ^7 and which 
describes the work rate on a stationary mass . The dis- 
placement work rate is identically equal to zero; this 
follows as a direct consequence of basic assumptions 3 
and 8: the absence of a potential field in which the 

subsystems move. 

Complete Determination of the System 

When the system is viewed as a whole, it is clear that the 
only input work is the work done by the adiabatic plane on 
the whole gas, the only output heat transfer is conduction 
at the head, and the difference must be the rate of change 


of internal energy of all the gas laminae. Expressed 
in equation form, it must be true that: 


_q wall (t) 



total (t) ” V total (t) 


dP(t) 

dt 


( 70 ) 


The system is correctly determined, then, if the N+l 
guiding equations contain Eqn. 70 as a corollary. 


The N first law equations for the N laminae can be 
written in the form of Eqn. 50: 


Qi (t) 
x in 



q wall (fc) 



(t) 



(t) 


= me 


* me 


= me 


dT 1 (t) 


mRT i (t) op(t) 

dt 


L P(t) dt 

dT n (t)‘ 


mRT n (t) dp(t) 

dt 


P(t) dt 


' dT N (t) l 

mRT N (t) 

dP(t) 

dt ^ 

j 

P(t) 

dt 


(71 

A,B,C) 


If we add these N equations, noting that the conductive 
flux into the nth element is the conductive flux out of 
the (n+l)st , the result is a gross cancellation on the 
left side, yielding: 


" q wall (fc ) 


me 

P 


2 


dT n (t) 

dt 


mR dP(t) 
PTtT‘~dt 


J. 

> 

n*l 


(72) 


The first term on the right side is really: 


1 


82 


mc„ 


N 

r 


n»l 


dT n (t) 

dt 


Y^' dH r. (t> 

Z , dt 

n=l 


dt 


(73) 


The last term on the right side is simplified with the help 
of the other basic governing equation, Eqn. 46: 


V d (t) ■ 2Z V n (t) ■ 

n=l 


N 



n=l 


mR 

pCt7 



(t) 


mR 

FTt T > 


T n (t >. 


(74) 


Then the last term in Eqn. 72 is: 


mR dP( t ) 
P(tT dt 


N 

(t> - 

n=l n 


v dP(t) 
d dt 


(75) 


Having re-expressed the terms of Eqn. 72, we collect them: 


q wall^ 


dH total (t > 

dt 


v dp ( fc ) 
a dt 


(76) 


This is the same as Eqn. 70 j thus, the N+l chosen equa- 
tions correctly determine the system. 


SOLUTION OP THE FINITE DIFFERENCE EQUATIONS 

The unknowns in the defining equations, Eqns. 46, 61, 64, 
and the N-2 versions of 60, are the future temperature 
array T n (j+i), n ■ 1,2,3»...N and the future chamber pres- 
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sure P(j+1). The method of solution chosen is as 
follows : 

1) Select the Initial temperature array T n (J=0) 
and the initial pressure P(j=0) from experience. 

2) Guess the future pressure at time 1 from the 
polytropic relationship: 

’V,(j=0)1 n 

P(j=i) “ P(j=0) v ^ QVO ( 77 ) 

where n has some reasonable value, from 1.35 to 
1.40 (in the program, n = 1.4 was chosen arbi- 
trarily) . 

3) Each of the N energy rate finite difference 

equations has two unknowns: the pressure P(j+1) 

and the temperature T n (j+1). The pressure having 
been guessed, each of the energy equations can be 
solved for the temperature. This yields the 

' ' ' • . . _ i .V , 

temperature array T n (j=l), n = 1,2,...N, for the 
guessed pressure at that time. 

4) The volume constraint equation, Eqn. 74, is 
written in the form: 

N 

P better^ -1 ^ 
guess 



and the calculated temperature array is inserted. 
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producing the P better EUCSS consistent with 
the constraints. This P(j=l) is compared to 
the guessed pressure from step 2. If it agrees 
within .002 psl (generally it agrees within about 
0.5 psi on the first iteration), it was a close 
enough guess . If not , this new pressure Is used 
in step 3 and the iteration begins anew. Conver- 
gence generally requires from three to six Itera- 
tions . 

5) In this fashion the temperature array (or profile) 
and gas pressure are calculated around the cycle 
from the "experience-selected" initial conditions. 
The temperature array and pressure value at the 
cycle end (-180 crankangle degrees, or BDC) are 
compared to the initial conditions. If the 
pressure has converged to within 0.1 psi, the 
cycle is deemed converged. If the pressure has 
not yet converged, the final conditions are used 
as initial conditions in step 1 and we begin 
again. With judicious selection of initial con- 
ditions, about five or six cycles are needed to 
achieve complete convergence* 

The computer program as written required about 27 seconds 
for compilation, and about two minutes, 12 seconds per 
cycle on the CDC 1604 computer, when N is set at about 70 
elements. Thus, the program was relatively expensive to 
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run. More sophisticated programming could reduce the 

time somewhat, but the very large number of Iterations 

required dictates that the execution time will be fairly 
long. 

CONVERGENCE AMD STABILITY 

With this program as with all finite difference treatments 
of differential equations, the matters of convergence and 
stability are of real concern. Two basic parameters can be 
varied independently to change the rate and/or accuracy of 
computation: the width of the spatial element, Ax n , and 

the time increment between steps. At. In the present case, 
the time increment is related to system frequency, rpm, and 
the time division of the individual cycle (here taken as one 
crankangle degree, with 360 degrees per cycle). The spatial 
width Ax n of the elements varies with the element chosen 
and also with time. 


The stability criterion for the Diffusion Equation is gen 
erally taken to be: 


k At 
gas 

pCp( Ax)‘ 


0.5 


( 79 ) 


However, p is varying from place to place and with time, k 
and c p are temperature-dependent, and the Ax of the elements 
are different and time-varying. In the face of these fac- 
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tors, a more pragmatic approach to stability was adopted. 

The time increment can be written (having fixed 360 as the 
number of increments per cycle) as 

At = C-^/rpm, 

where is some constant, and the Ax increment can be 
roughly taken to be: 

x = C 2 /N, 

where C 2 is some number and N is the total number of 
laminae. Then in an approximate fashion, the stability 
requirement can be expressed as : 

N 2 

— < some number C, @0) 


where C is determined by trial and error. The table below 
gives the results of several trials: 


STABILITY CORRELATION TABLE 


Run Number 

RPM 

N 

Calculated C 

Solution Condition 

1124 

400 

42 

4 .410 

Unstable 

1187 

400 

40 

4.000 

Unstable 

621 

600 

71 

8.402 

Unstable 

993 

600 

35 

2.042 

Stable 

1064 

800 

50 

3.125 

Unstable 

890 

800 

42 

2.205 

Stable 

629 

4000 

150 

5.625 

Unstable 

544 

4000 

71 

1.260 

Stable 

784 

9000 

71 

0.560 

Stable 
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As N is decreased below about 30, the Ax "screen” gets 
rather coarse, and the program results become* a less faith- 
ful solution to the differential equations. 


MODEL RESULTS AND CONCLUSIONS 

As this model differs from the experimental apparatus in 
a number of ways, care must be exercised in making direct 
comparisons between model and experimental dependent varia- 
ble arrays. The one-dimensional assumption neglects heat 
transfer to the cylinder sleeve. The presence of viscous 
energy dissipation and very small spatial pressure gradients 
is ignored. In the experimental setup, the energy dissipated 
at the piston rings and piston surface in rubbing on the 
non-lubricated cylinder walls was thought to contribute in 
determining the general thermal level of the working fluid. 
This energy dissipation is not included in the present 
model. The mass of the model gas is constant; in the exper- 
imental setup, some mass was lost through blowby during 
the high pressure portion of the cycle and regained through 
the makeup ports in the low pressure portion of the cycle. 
This fluid flow is analyzed in the system simulation pro- 
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gram (SIMUL) , not presented here. A final impor- 
tant difference is the lack of makeup air in the model, 
which contributed a fresh "charge" of 140°F air to each 
experimental cycle. This had two effects: first, a cer- 

tain amount of energy was added to the system with the 
makeup air; second, enough mass was added to keep the 
pressure near BDC at a level near atmospheric, from which 
compression could proceed. 

Table 1 on the following two pages lists the independent 
and dependent variables for the model as a function of 
frequency; a few parallel experimental variables are in- 
cluded for purposes of comparison. Peak crankangle times 
listed in the table are interpolated between the computed 
one degree steps, and because the solution has a step, or 
segmented, nature, they should not be thought of as being 
accurate to the nearest 0.1 crankangle degree. 

The method of calculation of the bulk-averaged temperature 
is described later. The polytropic coefficients based on 
the bulk averaged temperature and the volume (I .365 at 
800 rpm, 1.370 at 4000 rpm) are in reasonable agreement 
with the experimental coefficients (1.359 at 915 rpm, 

1.374 at 3452 rpm). 

The information (temperature swings and ratios of max/ 
n*in) listed for elements 1 and 2 should be approached with 
caution, since the masses and volumes of the elements de- 


Table l 

CHANGE OF VARIABLES WITH FREQUENCY - 
ADIABATIC PLANE MODEL AND EXPERIMENTAL DATA 


ADIABATIC PLANE MODEL 


RPM 

600 

800 

4000 

9000 

Frequency, Hertz 
N, no. of elements 
Chosen Wall Temp., °F 

10 

35 

200 

13.3 

42 

215 

66.7 

71 

276 

150 

71 

276 

Swines in Variables 





T (bulk averaged), °F 
T (first element), °F* 

T (second element), °F"" 
T (Nth element), °F 
Pressure, psi 
q'at wall, Btu/hr ft 2 

126.1 
293.0 
624.5 
130.2 
11. 9K 

530.0 

120.5 

284.7 
66 5.1 

133.8 
14 .OK 

665.0 

149.8 
343.3 
850.0 

152.9 
33. 9K 

216.9 

445.5 

891.1 

159.4 

52.6K 

Ratios of Max/Min 





T (bulk averaged) 
T (first element) 
T (second element) 
T (Nth element) 

1.218 

1.616 

2.505 

2.350 

1.202 

1.575 

2.500 

2.380 

1.234 

1.655 

2.442 

1.360 

1.936 

2.434 

Pressure 

24.81 

24.95 

25.16 

25.22 

Pressure at BDC, psia 

5.47 

5.59 

6.34 

6.58 

EXPERIMENTAL COUNTERPARTS. 



RPM 


915 

’ 3452 


Frequency, He r t z 
Measured Wall Temperature ° 

F 

15.27 

225 

57.5 

277 


Swings in Variables 





Gas Temperature (from wire) 

Gas Pressure, psi 

Wall Heat Flux (Btu/hr ft 2 ) 

, F° 

560 
365.4 
45. 3K 

630 
364.0 
94. 5K 

' 

Ratios of Max/Min 





Gas Temperature 
Gas Pressure 


2.113 

23.9 

2.188 

30.5 


Pressure at BDC. psia 


17.61 

12.04 



"element next to the wall 
""second element from the wall 


Table l (continued) 


90 


ADIABATIC PLANE MODEL 


RPM 

6 oo 

800 

4000 

9000 

Frequency, Hertz 

10 

13.3 

66.7 

150 

Peak Crankangle Times* 

T (element N),cr. degrees 

- 0.2 

- 0.1 

+ 0.2 

+ 0.2 

Pressure, » 

- 0.3 

- 0.2 

- 0.1 

- 0.0 

q at wall surface, <■ 

- 10.2 

- 10.3 

- 9.8 

- 9.6 

Polytropic Calculations 
n from T bulk- v total 


1.365 

• 1.370 


n from P, V tot<a 

1.392 

1.392 

1.382 

1.381 

Cycle averaged T^, °F 

134.5 

174.7 

372.3 

441.4 

Cyclic mean y(T N ) 

1.398 

1.396 

1.388 

1.386 

Ratio of n(P,V tot )/y 

0.982 

0.986 

0.993 

0.996 

EXPERIMENTAL COUNTERPARTS 



RPM 


915 

3452 


Frequency, Hertz 


15.27 

57.5 


Peak Crankangle Times 
Gas Temperature 


- 2.0 

+ 2.2 


Gas Pressure 


- 1.0 

- 1.0 


q at wall surface 


- 8.5 

' - 8.5 


Wall Surface Temperature 
Polytropic Calculations*** 


+ 67 . 5 "* 

- 2.4 


n fr0l “ W V total 


1.320 

1.335 


n from p 6 aa> Y total 


1.359 

1.374 



"Computer iterations are at one crankangle degree Inter- 
vals, and values listed in tenths of a degree represent 
interpolations . 

""This is the late "hump" which occurs only at low speeds. 

"""The experimental polytropic calculations are performed 
for a system whose mass is time-varying, and are to be 
interpreted with caution. 
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pend on N, the number of elements chosen. As this varies 
with frequency (because of stability requirements), these 
data are strictly interpretable only within a frequency 
(or rpm) column. Note that the temperature swing of the 
first element is generally about half that of the second 
element, and is on the order of one or two hundred degrees. 
If we would wish the model to precisely describe the finer 
details of the temperature distributions, we would have to 
increase N to several hundred, greatly increasing the com- 
puting time. Since temperature ratios on the order of 2.5 
represent isentropic behavior and a unity temperature ratio 
isothermal behavior, it can be seen from the listed temper- 
ature ratios how we proceed from isentropic behavior near 
the Nth element to nearly isothermal at the first. 

Pig. 25 shows the histories of the adiabatic plane temper- 
ature (actually the temperature T^, half the thickness of . 
the Nth lamina from the adiabatic plane and never more than 
0.2 F° from the plane temperature) with frequency as the 
parameter. This temperature peaks consistently within ±0.2 
crankangle degrees of TDC . The swing and level increased 
with increasing frequency, and the temperature ratio de- 
creased very slightly. The shapes of these curves and 

tne crankangle locations of their maxima compare well with 

■) 

most of the experimental observations. The polytropic 
coefficient calculated from the T^ ratio and the overall 
volume ratio decreased from 1.3922 at 10 Hz to I. 38 O 8 at 
150 Hz. 


temperature 
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Pig. 26 gives several gas pressure histories as pre- 
dicted by the model. The pressure was closely symmetric 
with crankangle in the compression and expansion strokes, 
and peaked within a few tenths of a crankangle degree of 
TDC at all speeds. The pressure swings and levels dropped 
as speed decreased below a few thousand rpm. As the lowest 
pressure in the cycle, (i.e., at BDC ) dropped from 6.58 psia 
at 150 Hz to 5.47 psia at 10 Hz, the cycle pressure ratio 
remained fairly constant at from 24.8 to 25.2. The poly- 
tropic coefficient calculated from the pressure ratio and 
the overall volume ratio increased from 1.373 at 10 Hz to 
1.3808 at 150 Hz. The pressure ratios compared well with 
the experimental and SIMUL values . The pressure swings 
were reduced because in the model the BDC pressures were 
much lower than in the experiment. The minimum pressures 
in the model (and hence the pressure and temperature 
swings) could easily be elevated by increasing, the mass 
of the elements, which is established in the initial part 
of the program. In the experimental apparatus, this was 
done in the introduction of makeup air. 

It is interesting to note the polytropic calculations for 
the Adiabatic Plane Model on the top of the second page of 
Table 1 . The cycle averaged adiabatic plane temperature 
rises from 134. 5°P to 44l.5°P from 10 Hz to 150 Hz. The 
ratio of specific heats, y» evaluated at the cyclic mean 
temperature, decreases from 1.398 at 10 Hz to I .386 at 150 



CRANKANGLE 
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Hz. The "degree of isentropicity " is best Judged from 
the ratio of the cycle polytropic coefficient based on P 

and V total* t0 The °y cles appear to approach the isen- 
tropic standard as frequency increases; this agrees with 
intuitive reasoning, in that heat transfer should be of 
less importance at high speeds, when the time interval for 
transferring heat per cycle is small. 

Fig. 27 presents the wall surface heat flux histories as 
a function of frequency . The histories reached their peaks 
at -10.0 ±0.iJ crankangle degrees, in close agreement with 
the high frequency experimental data. The history shapes 
were very close to those of the high frequency ( >25 Hz) 
experimental heat flux histories (see Fig. 13 for experi- 
mental plots), although the amplitudes were slightly 
smaller. Since the surface temperature maxima were found 
in experimental and analytical studies to lag the corres- 
ponding surface heat flux maxima by from 9 to 11 crankangle 
degrees for this shape of history, the transient surface 
temperature which could be calculated from these model 
heat flux histories would peak at ±1.6 degrees from! TDC. 
This correlates very closely with the high speed experi- 
mental observations. 

Figs. 28 and 29 give the temperature profiles in the gas 
between the head and the adiabatic plane. Observe that a 
true boundary layer thickness would be somewhat difficult 
to evaluate from these data, although the steepest temper— 


HEAD SURFACE HEAT FLUX 



CRANKANGLE 






DISTANCE INTO GAS FROM WALL SURFACE, inches 

Pig. 28 
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ature gradients seem to be within 2mm of the wall at BDC 
and 1mm at TDC. The traces for +20 and +40 crankangle 
degrees on the expansion stroke (Pig. 29 ) show clearly 
the "negative h" phenomenon: with a positive temperature 

potential from gas to wall, the heat flux is from wall to 
gas. The +20° trace corresponds to a zero or slightly 
negative value for h; a +60° trace (large heat transfer, 
nearly zero AT) would correspond to nearly infinite h. 

Fig. 30 shows how the profile at 30° before TDC (chosen 
only because of its distinctive shape) changes with 
driving frequency. It is clear that the boundary layer 
thickness is much larger at low frequencies than at high 
frequencies. 

Fig. 31 shows the behavior of the heat transfer coeffi- 
cient calculated in the program from the ratio of the wall 
heat flux and the adiabatic plane to wall temperature dif- 
ference. This coefficient is compared to information pub- 
lished by Overbye (33) which was generated in a motored 
CFR engine with four piston strokes per data cycle and 
which required a calculated gas temperature. It is inter- 
esting to note that in spite of the difference in engine 
cycles and in the orders of magnitude of the heat transfer 
coefficients, they pass to infinity (i.e., the gas - wall 
temperature difference passes to zero) at nearly the same 
times and have quite similar shapes. The vertical arrows 
on the abscissa indicate the infinity asymptotes. The dif- 






TRANSFER COEFFICIE 


HEAT TRANSFER COEFFICIENT (FOR OVERBYE, Btu/hr ft ) 
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CRANKANGLE 



102 

ference In the sign of the histories after +100 crankangle 
degrees is due to the fact that in the Overbye data, the 
heat transfer was still from gas to wall at this stage, 
whereas the model predicts small negative heat transfer. 

It is possible to calculate local fluid velocities from 
the rates of change of location of the laminae in the gas. 
When this was done, it was found that the velocity profiles 
from piston to wall were very close to linear for most 
portions of the calculated cycle. The nonlinearities 
occurred generally near the wall and were most obvious near 
TDC and BDC: where the adiabatic plane velocity changed 

sign precisely at -180° and 0°, the velocity at the center 
of the first lamina switched sign at -172° and +2° at all 
frequencies. In the larger portions of the compression 
and expansion strokes, however, the nonlinearity amounted • 
to a few per cent of the velocity value calculated from 
the linear assumption. 

A bulk-averaged gas temperature could be calculated at 
any instant in time by integrating the temperature profile 
in the gas with respect to distance. An easier, but equi- 
valent, method is to use the perfect gas law: 

P V 

m „ total 
bulk mR 

A comparison of the bulk gas temperature to the adiabatic 
plane temperature at 4000 rpm is presented in Fig. 32 . The 
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bulk temperature' is seen to be consistently below the 
adiabatic plane temperature. 

Since a heat transfer coefficient is generally applied 
to a bulk gas - wall temperature difference, a coefficient 
history can be calculated for the cycle from the relation- 
ship 

h adiabatic plane model a q wall^ T bulk “ T wall^’ 

For purposes of comparison, the Eichelberg coefficient can 
also be calculated from the driving frequency and the pres- 
sure and bulk gas temperature histories. These coeffi- 
cients are compared in Fig. 33 . 

A comparison of the several approaches to wall heat trans- 
fer is presented in Fig. 3 ^. The adiabatic plane model 
heat transfer is seen to agree quite well with the experi- 
mental heat flux with regard to phasing and shape. The 
Eichelberg heat transfer maximum is 12 degrees late and 
the negative heat transfer portion of the cycle is poorly 
simulated. 

Finally, Fig. 35 shows the poly tropic coefficient histories 
for the Nth volume (next to the adiabatic plane). These 
coefficients were calculated at time J from the data at 
time j-1 with the relationship: 


n N (J) “ VCl-(log[T N (j)/T N (j-l)3/logCP(J)/P(j-l)])]. 
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Fig. 33 
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(for Eichelberg and Model Traces only) 
HEAT TRANSFER COEFFICIENT, Btu/hr °F ft 2 




COEFFICIENT HISTORY 
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Fig. 35 
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The coefficient for the Nth element is seen to be close 
to isentropic throughout the cycle for all speeds except 
near BDC and TDC (i.e., ±180°, 0°). The anomalous be- 
havior at these points is probably due to the effect of 
the step nature of the calculation at the points of zero 
rate of change of the variables. 

It was hoped that a generalized correlation could be de- 
veloped for the prediction of accurate heat transfer his- 
tories from relatively simple engineering data: mean wall 

surface temperature, bulk gas history, and Eichelberg 
coefficient. Much effort was expended toward this goal. 
However, the variable histories did not have an easily ex- 
pressed frequency dependence. The surface heat flux maxi- 
mum and minimum occurred at the same crankangle position 
over the frequency range, but the crankangle times at which 
the difference (T^^ - T wa -^) passed through zero changed 
with frequency. As a consequence, the ±» asymptotes of the 
film coefficient were frequency-sensitive; the film coeffi- 
cient could not be expressed as a simple function of fre- 
quency, wall temperature, and instantaneous bulk properties. 

In conclusion, this model for the working fluid is pre- 
ferred over the so-called secondary models, in that it 
proceeds from basic thermodynamic principles, is guided by 
a few simple assumptions, and contains little or no empir- 
ical data. It is interesting to note, however, that many 
of the empirical estimations used in earlier models agree 
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reasonably well with the calculated behavior of this model', 
and the output data (primarily the wall surface heat flux 
history shapes and phasing) of the Simulation Transient 
Model, the Variable h Resistive Model, and this Adiabatic 
Plane Model are in fairly good agreement. 

The Adiabatic Plane Model output is in very close agree- 
ment with the high frequency experimental data with res- 
pect to phasing and the dependent variables' history shapes, 
although the variable amplitudes are slightly smaller, proD- 
ably because of the influence of friction and makeup air 
in the experimental system. The program for the model, 
although presently rather expensive in terms of computation 
time, is very flexible. An interested investigator with 
sufficient time could modify it to: 

1) study the effect of working fluid properties on 
heat transfer; 

2) include time-dependent energy generation in the 
different laminae to simulate combustion or dis- 
sociation; 

3) include the effects of viscous dissipation; 

4) include time-dependent terms for gas - cylinder 
wall energy exchange and piston friction energy 
dissipation; and 

5) compute the head surface temperature history from 
the surface heat flux history (as Fig. 27 , calcu- 
lated in program with constant surface temperature) 
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insert time-dependent wall temperature as new 
boundary condition and rerun program. Comparison 
of wall heat flux histories would provide a mea- 
sure of "feedback" from wall surface to gas. 


Ill 

IX. DISCUSSION OP RESULTS AND CONCLUSIONS 

Experimental bulk gas temperature histories are given in 
Appendix A. As the driving frequency increases from 10 to 
60 Hz, the crankangle location of the gas temperature maxi- 
mum shifts from -2.5 to -1.5 crankangle degrees before top 
dead center, probably because of a phase lag introduced by 
the resistance wire (see Section VII. A.). The gas tempera- 
ture swing increases from about 560°F to about 630°P. The 
histories remain approximately symmetric about the TDC or- 
dinate line. 

As the driving frequency increased from 10 Hz to nearly 
60 Hz, the swing of the experimentally measured gas tempera- 
ture increased from 560 to 630°P. By comparison, the 
Adiabatic Plane Model amplitudes increased from 625 to 
850°P from 10 to 67 Hz. This amplitude difference can be 
attributed partially to an unknown and uncalculable attenua- 
tion of the resistance wire transducer due to thermal 
inertia and thermal resistance in the wire’s own boundary 
layer, and partially to a slight mismatch of the Adiabatic 
Plane Model to the actual physical system (including, as 
previously noted, viscous energy dissipation and makeup and 
blowby air). None of these factors are very well understood 
quantitatively; however, small inaccuracies in the gas 
temperature amplitude have little effect on the qualitative 
side of the observations. The great difficulties in ob- 



112 


tainlng even a poor local gas temperature history were 
noted In a previous section. 

Experimental bulk gas temperature ratios (maximum/minimum) 
ranged from 2.113 at 15.25 Hz to 2.188 at 57.53 Hz, Adia- 
batic Plane Model gas temperature ratios (bulk averaged gas 
temperature) were 2.350 at 13.3 Hz and 2.380 at 66.7 Hz; • 
these ratios are somewhat higher than the experimental 
ratios because the model minimum temperature was somewhat 
low (because the mass in the tabulated example was guessed 
lower than the experimental mass— hence the compression 
process of the equilibrium cycle began at a low gas tem- 
perature and pressure). 

Gas pressure histories are also found in Appendix A. The 
history shape, amplitude, and location of maximum are 
essentially independent of driving frequency: the amplitude 

remains at about 360 or 365 psi, and the maximum occurs at 
-2.5 crankangle degrees before TDC over the frequency range 
10 to 60 Hz. 

The gas pressure histories predicted by the Adiabatic Plane 
Model are almost identical to their experimental counter- 
parts; the history shapes are much alike and the maxima 
occur at a few tenths of a crankangle degree before TDC. 

The amplitudes predicted by the model are low with respect 
to the experimental amplitudes because the mass assumed in 
the tabulated example is low. A later computer run of the 
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model with a higher system mass (or, equivalently, a more 
reasonable gas pressure at BDC of the equilibrium cycle) pre- 
dicted gas pressure amplitudes close to the experimental 
values. 

The Adiabatic Plane Model pressure ratios (maximum/minimum) 
were in the range of 24 to 25 on all runs. The experimental 
ratios ranged from 23.9 at 15.25 Hz to 30.5 at 57.53 Hz. 

The latter were subject to perhaps 103» scatter, however, 
because of resolution problems: an error of 1 % in the min- 

imum pressure (say 3.7 psi at 14 psi) causes a large change 

in the pressure ratio (from 26.9 to 33.4). Pressure ratio 

> 

is not, therefore, the best of parameters by which to des- 
cribe these pressure histories; pressure swing might be 
better. 

The experimental surface temperature histories can be 
found in Appendix A. Fig. 13 is a summary of surface tem- 
perature and surface heat flux histories, and with Fig. A1 
for a lower frequency constitutes the best illustration for 
the following discussion. The data of Fig. A1 were not 
included in the summary figure, because gas properties were 
‘not measured during this run; the runs of Fig. 13 represent 
complete data sets. Fig, 14 shows the locus of the surface 
temperature maxima with driving frequency, and will also be 
helpful. 

At low frequencies, the predominate characteristic of the 
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experimental surface temperature history is seen in essence 
in Pig. Al: the history is relatively smooth and nearly 

sinusoidal, peaking at +55 crankangle degrees at 9.65 Hz. 

As driving frequency decreases, the so-called "hump 1 ' shifts 
closer to the gas temperature maximum (near TDC, or zero 
crankangle degrees); no experimental recordings were made 
at less than 9*65 Hz, but at frequencies less than, say, 

1 Hz, the gas temperature and wall surface temperature his- 
tories should be nearly in phase. As the frequency in- 
creases to 15 or 16 Hz, the peak shifts later in time (to 
+66 degrees at 15.25 Hz, for example). Above this fre- 
quency, the late "hump" degenerates in amplitude and a 
second peak begins to emerge at TDC (refer to Fig. 14). 

At frequencies greater than about 25 Hz the "central" peak 
predominates. Above 33 Hz, the late peak is hardly dis- 
cernible and the surface temperature histories are somewhat 
peaked and nearly symmetric about the TDC ordinate line. 

These trends are reflected in the experimental heat flux 
histories of Pig. 13 (with Pig. Al for 9.65 Hz). At very 
low frequencies the heat flux histories are also roughly 
sinusoidal, leading the surface temperature by some 20 or 
25 crankangle degrees. As the frequency increases the two 
heat flux "humps" separate, the late hump growing smaller, 
the central peak growing and leading the surface tempera- 
ture maximum by about 17 degrees. At 57.53 Hz, the late 
hump in heat flux is quite small. The high frequency ex— 
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perimental heat flux histories are very much like those 
predicted by the Adiabatic Plane Model in Fig. 27. 


This is not the first occasion on which this type of be- 
havior in surface temperature or surface heat flux has 
been observed. T. LeFeuvre (80), working at the Univer- 
sity of Wisconsin, used a surface thermocouple to measure 
instantaneous head surface temperatures in a diesel engine 
under both motored and fired conditions. The single— cylin- 
der research diesel engine had a compression ratio of 15:1 
and a swirl ratio of about 2 (i.e., the working fluid 
swirls, or rotates, twice inside the compression chamber 
during one revolution of the crankshaft. Since this was a 
four-stroke cycle device, two revolutions of the crankshaft 
represent one cycle of the system) . Under certain condi- 
tions LeFeuvre obtained surface temperature histories with 
two maxima remarkably' like those previously presented (at 1 
about 16 tb l ‘20' Hz). “Fig. 3 6 is from ah oscillograph re- 


cord of LeFduvre’s 7 observations. “ Most of the "hash," or 
high frequency- Slghal, 1 on the temperature trace is probably 
aerodynamic “in L Oirigin 1 ,' since it was closely 1 reproducible* 
from cycle to 6jrcle (particularly in thb motored case), 
and probabl^ “represehted L Regular patterns of' air motion. 
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A.Chu ( 8l) reported unpublished work done at Cummins 
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experimental gas pressure histories were used to compute 
bulk gas temperature histories; the Annand heat transfer 
coefficient was used to represent the thermal resistance, 
but other details of the model used for the boundary layer 
were not specified. At 1500 rpm (25 Hz) the computed 
head heat flux history had a double peak, one at about zero 
crankangle degrees, the other quite late in the cycle. By 
2100 rpm (35 Hz) the history had resolved to the expected 
peak at about TDC. No further information was available 
on these observations. Chu also reported encountering an 
experimentally measured head temperature history showing 
the double humping phenomenon. 

It is difficult, if not impossible, to construct a single 
theoretical model predicting both the high frequency and 
the low frequency experimental behavior of the thermal 
boundary layer. Some modelling can be done, however, if 
the regions are considered separately. 

In the experimental observations at frequencies greater 
than about 25 Hz (see Fig. 14), the central peak is the 
strong feature of the surface temperature history; peaking 
at about zero crankangle degrees (i.e., the top dead center 
position of the piston), the surface temperature lags the 
gas temperature by only two or three degrees. The surface 
heat transfer histories associated with these surface tem- 
perature histories led the surface temperature by about six 
crankangle degrees and had the characteristic shape described 
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in Pig. 14. The Adiabatic Plane Model did a good job 
in predicting the high frequency heat transfer in that the 
shape and the phasing of the wall surface flux histories were 
quite close to their experimental counterparts. 

The low frequency surface temperature observations are more 
troublesome to understand. The RC model for the gas-wall 
interface, which is explained in detail in the author's thesis 
assumes that the thermal boundary layer is relatively thick 
(four millimeters on the average) and that pressure work 
exerts no influence in the boundary layer. This model will 
account for the low frequency (up to about 20 Hz) surface 
temperature and heat flux histories. However, as this model 
does not represent a smooth transition from a theoretical 
treatment of the energy equation to a mathematical represen- 
tation of the energy fluxes in the boundary layer, and is a 
rather after-the-fact construction, it is not presented in 
this report. 

The significance of the transition region (i.e., 18 to 30 Hz), 
in which the surface temperature history has two maxima, is 
unknown. It suffices to say that other investigators have 
observed this phenomenon and have offered no explanation 
for it. 

Finally, at this point we can summarize the important 
thoughts scattered in various places in this text on the- 
relationship of surface heat transfer response to the gas 
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pressure and temperature driving functions. The lag or 
lead of surface heat transfer at the surface of a fuel drop- 
let with respect to variations of pressure and/or tempera- 
ture with time is a crucial factor in determining whether 
increased vaporization of the droplet will support or oppose 
the amplitude of the variations. The magnitude of this lag 
or lead is dependent upon, among other things, the proper- 
ties of the gas and the wall, and the harmonic 'distribution 
of the pressure and temperature oscillations. 

We can focus our thoughts by "constructing” a gas in the 
one-dimensional space between a plane wall surface and a 
bulk gas "driver" which is undergoing large periodic excur- 
sions in pressure and temperature. If we begin by hypothe- 
sizing a gas with no heat capacity and an invariant pressure 
but only a thermal conductivity (which can be thought of as 
indicating only a thermal resistance between wall and driver 

gas), the driving potential for heat transfer fT (t) - 

gas 

T wall^^ and the conductiv ity of the gas (or the reciprocal 
resistance) will maximize at approximately the same time as 
the "driver" temperature. As a consequence, the surface 
heat transfer will be almost exactly in phase with the tem- 
perature and pressure of the driving gas. The obvious ana- 
log is a resistive electrical circuit. The comparable model 
in this work is the Variable h Resistive Model, and as ex- 
pected, the heat transfer peaks right in phase with the 
gas temperature (and with the gas— to— wall temperature dif— 



ference) at about -1 crankangle degree. 
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If, still keeping the pressure constant, we add to the gas 
a finite heat capacity, the heat flux and temperature varia- 
tions of the surface are delayed with respect to the temper- 
ature variations of the bulk gas driver: when the driving 

temperature is rising, the capacity causes the interface gas 
to act as an energy absorber, decreasing the energy flux to 
the surface. When the temperature falls, the interface gas 
releases the energy. The net result is a lag, or delay in 
the peaking of the surface temperature and heat flux with 
respect to the driving gas temperature. The RC Model in 
this work examined this phenomenon by varying the mean value 
and the amplitude of the heat capacity of the interface gas 
while keeping the mean value and the amplitude of the ther- 
mal conductivity the same. As an example, in a computer run 
at 915 rpm (15.25 Hz) with a Nu^ of 15 and a mean boundary 
layer thickness of 4.6 mm, the heat flux peaked at 40 degrees 
after TDC (or 42 degrees after the bulk gas temperature). 

When the heat capacity was decreased to zero, the heat flux 
shifted and became exactly in phase with the bulk gas tem- 
perature again. 

When the pressure of the interface gas is forced to vary 
with time (and essentially in phase with the bulk gas tem- 
perature), the pressure work in the interface gas leads the 
bulk gas temperature and pressure variations. The phase 
of the heat transfer with respect to the bulk gas tempera— 
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ture variations is then determined by the relative magni- 
tudes of the pressure work and the capacitive energy storage 
In the experimental measurements with air, the heat flux 
peaked at about 8,5 crankangle degrees before TDC (about six 
degrees before the bulk gas temperature and pressure). In 
the Simulation Transient Model, the heat flux peaked at 
about 6 to 10 degrees before TDC. In the Adiabatic Plane 
representation, the heat transfer led the bulk- gas tempera- 
ture and pressure by 9 or 10 crankangle degrees. These 
phase relationships are illustrated in the sketch below, 
which is not to scale. 



Thus, the '‘leading" effect of pressure work slightly more 
than offsets the "lagging" effect of gas capacity. It is 
apparent that pressure work is sin important factor in de- 
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termining the phase relationship of the heat transfer to 
the "bulk gas" oscillations. 

In summary, the following statements can be made from the 
experimental and theoretical observations: 

1) Above about 30 Hz (to about 60 Hz experimentally 
and to at least 150 Hz analytically) the wall 
surface temperature lags the bulk gas ‘temperature 
by only 2 or 3 crankangle degrees (.03 to .05 
radians) and the gas pressure by 1 or 2 degrees 
(.02 to .03 radians). 

2) The surface heat transfer rate leads the surface 
temperature by about 10 crankangle degrees. 

3) The surface temperature history shape or ampli- 
tude do not change appreciably above 30 Hz. 

4) The amplitude of the heat flux history increases 
regularly with frequency, from 52,400 Btu/hr ft 2 
at 25.21 Hz to 93,900 Btu/hr ft 2 at 57.53 Hz. 

5) For frequencies above about 30 Hz, a simple re- 
sistance model of the boundary layer (with a 
time-varying resistance based on a "reasonable" 
internal heat transfer coefficient) gives a 
satisfactory first order prediction of the thermal 
response of the wall to a gas-side temperature 
disturbance. The Adiabatic Plane Model, which 
does not require such an assumed coefficient, does 
an even better job. 
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6) For frequencies below about 20 Hz, the RC Model 
can predict the phase lag of the surface tempera- 
ture with respect to the gas temperature. 

7) The heat transfer coefficient h, calculated from 

the ratio of q _ and (T „ - T ni ), assumes 
^surface gas wall * 

negative, zero, and infinite values, and as such 
is of questionable utility in accurately des- 
cribing the instantaneous surface heat flux. 

8) Among the models developed, the Adiabatic Plane 
Model, based almost entirely on fundamental 


thermodynamics, provides the best qualitative 
model for unsteady heat transfer between a plane 
surface of a solid and a gas experiencing pres- 
sure and temperature fluctuations. 
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APPENDIX A 

GRAPHS OF EXPERIMENTAL DATA 
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III. PRESSURE RECORDING CIRCUITRY AND NOTES 

The piezoelectric Kistler 601K pressure pickup generated a 
charge of 1.078 picocoulombs per psi of applied pressure. 

This change in charge was converted to a voltage signal by 
a Kistler #1713 charge amplifier in ‘short time constant 
operation, and further amplified by an Astrodata direct 
coupled amplifier. The signal was displayed oh one channel 
of the four-channel Tektronix 565 oscilloscope. Provision 
was made for injecting a 1000 Hz square wave of 1.0 volt 
amplitude into a calibration input jack on the Kistler 
charge amplifier. By this means a simulated 100 psi square 
wave could be superimposed upon the pressure trace to aid 
in pressure signal calibration. 

The 601H transducer produced smooth, trouble-free signals 
through all test phases. Very little cycle-to-cycle varia- 
tion was noted in the high-pressure (say -90° < crankangle 
< +90°) half of the cycle. In the low pressure (up to per- 
haps 80 degrees to either side of bottom dead center) por- 
tions of the cycle, some cyclic variation was noted, perhaps 
inside an "envelope" of about 20 psi. The 601H transducer 
does introduce an erroneous pressure indication when the 
transducer surface is stressed by cyclic temperature 
changes. With a proportionality constant of about 3 psi/F°», 

"Unpublished data of J. Alyea, Department of Mechanical 
Engineering, University of Wisconsin 



154 


a surface temperature swing of 2 F° will Induce a false 
pressure reduction of 6 psi, a 1 . 1 % error compared to the 
3^0 psi pressure swing. Having no precise information of 
the phasing of this small error, it was ignored throughout 
the tests. 

Fig. 4 is a block diagram of the pressure signal circuitry. 
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IV. SURFACE TEMPERATURE CIRCUITRY AND NOTES 

There are two popular varieties of surface thermocouples 
in the United States: those manufactured by Nanmac Cor- 

poration*, of Indian Head, Maryland, and those available 
from Heat Technology Laboratory*, of Huntsville, Maryland. 
The Nanmac thermocouples have a wider range of application 
than the second variety. The active surface can be eroded 
by combusting gases or rubbing friction. For a variety of 
reasons, however, the product of the second manufacturer 
was chosen. 

The surface thermocouples used in this study were the same 
type as were used by Pless ( 3 * 0 , Bennethum ( 35 ), and Over- 
bye (36) in recording surface temperatures in spark- and 
compression-ignition engine combustion chambers. Developed 
by Bendersky ( 37 ) in the early 1950's, the instrument is 
little changed from the original in design. Iron-constan- 
tan and iron-nickel metal combinations were obtained and 
mounted in a cast iron head for all major tests. For some 
preliminary studies in which a variety of head materials 
were used, an aluminum-nickel surface thermocouple was the 
transducer. 

The thermocouple construction was as follows: a nickel 

(or constantan) wire has its surface oxidized; this nickel 

*See List of Manufacturers and Addresses 
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oxide is an electrical insulator. The iron is drawn 
down on the nickel wire by a wire-drawing die. After 
metallurgically grinding and polishing the end of the re- 
sulting concentric cylinders, a one-micron-thick layer of 
constantan is plated on the end. Thus the thermocouple 
junction is located at the washer-shaped constantan-iron 
interface. Originally, nickel was used as the plating 
material, but since the one-micron plating is at essen- 
tially the same instantaneous temperature throughout, the 
law of intermediate metals does not require the plating 
and wire materials to be common; a better bonding substance 
was chosen in the constantan. The cylinder of iron is 
finally given a UNC 2-56 thread. 

Fig. 5 shows a thermocouple as it was mounted in the head, 
and Fig. 6 gives, in block form, the circuit in which it 
was the signal source. The gas-metal interface temperature 
is seen in the voltage between ground (G) and the oxided 
nickel wire (Nl) with its cold junction, and is indicated 

by E i n terior* The head exterior temperature (a constant 
when the system is in equilibrium) is seen in the voltage 
between ground (G) and a nickel wire (N2) held between the 
head and an iron mounting nut, with its cold junction, the 
voltage indicated by the symbol E ex terior* 

Since the interior temperature signal must be displayed- at 
high gain for good resolution, and since this signal con- 
tains a large dc component (the swing of 2 F° superimposed 
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Fig. 5 Thermocouple Construction and Hookup 












160 


upon a 200 P° surface-to-ice bath difference), it is much 

more convenient to display the voltage corresponding to 

the temperature difference through the head (E., , 

interior 

E exterior^* An indic ation of equilibrium is given when 
this displayed voltage difference ceases to increase. The 
voltage (N1 - N2) is the same as voltage (II - 12) since 
the only difference between them is two iron-nickel junc- 
tions at 32°P connected in opposition; the former signal 
is the one displayed. Because even at equilibrium the 
surface temperature (with a swing of perhaps 2 P° per 
engine cycle) appears to "hunt" or drift up and down by 
a degree or two in a periodic fashion at a fraction (like 
one-tenth) of the system fundamental, photographic records 
are made with a blocking capacitor feeding only the ac 
signal to the oscilloscope. 

The voltage difference is amplified 3000 times by an Astro- 
data d.c. amplifier and displayed at 100 millivolts/centi- 
meter on one channel of the oscilloscope. Because one 
side of the transducer is electrically grounded (by its 
design), 60 Hz (and higher harmonic) noise causes great 
trouble in this type of work, and efforts must be made to 
minimize it. All connections must be shielded; cables must 
be coaxial and shields must be grounded as close to a cen- 
tral point as possible to avoid ground loops. D.c. ampli- 
fiers must have differential inputs and outputs, and input 
and output grounds should be independent from each other 
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and from the case. The oscilloscope should be operated 
in differential mode. A. S. Buchman ( 38 ) presents a good 
discussion of this kind of noise abatement technique. Data 
were taken when the local electromagnetic noise level was 
loWj with all nearby fluorescent lamps and unshielded trans- 
formers turned off. Batteries were used for all possible 
power supplies (as opposed to ripple-containing rectified 
power supplies). The very high frequency noise level 
(contributed primarily by the oscilloscope amplifiers) and 
the 60 Hz noise level were held below approximately 0.18 P° 
and 0.12 F° equivalent temperature signal, respectively. 

Some of the surface thermocouples used in this study were 
originally purchased in the late 1950's and were used by 
Overbye (33) in a fired spark-ignition engine. They were 
recalibrated and cleaned with denatured alcohol. At speeds 
in excess of 1800 rpm (30 Hz) some of these older thermo- 
couples occasionally gave an output which was clearly not 
a temperature response. This extraneous signal was clearly 
recognizable by its appearance: a high amplitude, high 

frequency harmonic signal superimposed in "packets" of 
10 or 20 crankangle degrees length upon the temperature 
signal. High speed drum camera films showed these extra- 
neous signals to be of 3000 Hz frequency. The portions of 
the cycle in which these signals appeared remained the 
same from cycle to cycle. No mechanism for this anomalous, 
signal is proposed, except perhaps vibration of a partially 



detached one-micron film. No traces containing this 
type of signal were used or analyzed. Surface thermo- 
couples sometimes degenerate to this type of signal be 
fore they break up when used in a fired engine . 
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V. GAS TEMPERATURE MEASUREMENT 


HISTORICAL SURVEY 

Probably the greatest problem in contemporary engine 
(or compressor) instrumentation is the accurate measure- 
ment of gas temperature at a point. An instrument must 
be small enough not to influence the flow about it mater- 
ially, rugged enough to withstand high gas velocity and 
temperature and pressure extremes, and sensitive enough 
to detect changes in gas temperature at relatively high 
frequency (the 50th harmonic at 3600 rpm is at 3000 Hz). 

The alternative to this set of requirements is a computa- 
tional method: assume constant mass, poly tropic processes 

at fixed polytropic coefficient, and a perfect gas whose 
properties are time functions only. A new analytical tool 
is the thermodynamic simulation program, made possible by 
the advent of large-core-storage, high-speed digital com- 
puters. A program of this sort is discussed later. How- 
ever, in general, a mass of experimental data is given 
more credence than an analytical conclusion in a flow system 
of this order of complexity. 

Historically, many investigators have measured temperatures 
in operating and motored engines. Some of the earliest 
methods suggested and tried involved the use of either re- 
sistance thermometers or thermocouples. Coker and Scoble 
(39) used thermocouples made of flattened platinum-alloy 
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wires 500 micro-inches thick (12. 7y). Burstall used 
resistance thermometers made of 38 to 64 micron platinum 
to obtain temperatures in an operating gas engine. Schultz 
(40) used 2 micron flattened platinum wires in measuring 
fluctuating gas temperatures. The slow response rate, oxi- 
dation of the wires, and early wire breakage left much to 
be desired. Thermocouples rugged enough to withstand en- 
gine turbulence proved to be too large for rapid tempera- 
ture response. 

The work of Semenov (41) represents one of the most suc- 
cessful attempts to measure .instantaneous local gas temp- 
eratures and velocities in a motored piston engine with 
small-wire instrumentation. This work was done in the 
Soviet Union prior to 1959 and was translated into English 
in 1963. Semenov's interest was primarily in gas turbulence; 
his three basic objectives were: 

1) to study the variation of turbulence parameters 
throughout the motored cycle, 

2) to study the nature of the gas flow during the 
intake and compression processes, and 

3) to study the effect of certain operating variables 
on the level and the frequency distribution of 
the turbulence. 

Measurements were made in a single cylinder CFR (variable 
compression ratio) engine with a cylindrical compression 
chamber (with large overhead intake and exhaust valves) at 
low speeds (from 10 to 20 Hz) and at moderate compression 
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ratios (4 to 9.5) . 

Although Semenov's major objective was velocity and tur- 
bulence level measurements, local gas temperatures were 
measured in the course of his work. A constant tempera- 
ture hot-wire anemometer (see next major section for des- 
cription of this type of transducer) was designed; it 
consisted of an 11 micron diameter tungsten filament for 
instantaneous gas velocity and turbulence measurement, and 
a long (about one inch) 11 micron diameter tungsten fila- 
ment heated by a "weak current" for temperature compensa- 
tion purposes. Since the constant temperature anemometer 
is sensitive to both velocity and temperature changes, the 
latter filament tries to follow the gas temperature; the 
temperature compensation in the anemometer is electronic 
rather than computational. In order to locate the compen- 
sation wire so that it sees a gas temperature history 
similar to that experienced by the velocity sensor, a gas 
temperature survey was made throughout the chamber around 
the cycle with an uncompensated 8 micron resistance ther- 
mometer at a 15 Hz driving frequency and 6:1 compression 
ratio. 

A number of interesting observations were made by Semenov 
which may find application in the present work. It must be 
emphasized, however, that because of the intake and exhaust 
processes, the gas velocity and turbulence levels in the 
Semenov investigation probably were considerably greater 
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than those encountered in this work. Velocity measure- 
ments were made between the exhaust and intake valves; the 
valves represented some h 8% of the head area. The obser- 
vations of Semenov include: 

1) Immediately after closing the intake valve on 
the compression stroke, the mean gas velocity 
dropped from about 8 meters/sec. to about 1.5 
meters/sec. in 45 crankangle degrees. The 
"directional" gas flow quickly ceased and only 
turbulent eddies remained. 

2) Sharp velocity gradients were found only within 
about 3 millimeters of the wall; this may be an 
estimate of the viscous boundary layer thickness. ' 

3) Measurements with the 8 micron resistance ther- 
mometer indicated that the thermal boundary 
layer was from 1.5 to 2 millimeters in thick- 
ness (at 15 Hz) and that the gas temperature 
beyond this penetration was spatially uniform 
(but time-varying). This tends to confirm a 
model incorporating a "bulk gas" and a thin 
boundary layer. 

In the 1920's a trend toward using radiation techniques 
started. Hershey and Paton (42) of the University of 
Illinois measured temperatures using the Spectral Line 
Reversal Method. Rassweiler and Withrow of General Motors 
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(^3) used a similar technique. More recently, the tem- 
peratures in the unburned gas of an engine were measured 
by other optical methods. The absorption of radiation by 
iodine vapor was used in one technique (M) by Beck and 
Chen, and Agnew (^5) described a unit which operated by 
comparing the relative amounts of infrared radiation 
emitted in two wavelength bands. 

Another method of measuring fluctuating gas temperature in 
a cylinder-piston-compression chamber is to measure the 
time interval it takes for a high frequency sonic pulse to 
travel the bore of the cylinder. The computed speed of 
sound is then related to the mean temperature of a perfect 
gas. This method has been developed by C. P. Taylor (^6) 
of Sloan Laboratory, Massachusetts Institute of Technology. 
Still another is the infrared null method, used at the 
University of Wisconsin by Burrows (^7) to measure end 
gas temperatures in a spark ignition engine and gas tern- 
peratures up to and during injection in compression-igni- 
tion engines. Radiation from a blackbody held at an eleva- 
ted, fixed, and known temperature is passed through filters 
and sapphire windows across the bore of a compression/com— 
bustion chamber and allowed to fall on a lead sulphide 
sensor. When the gas is at the blackbody temperature, the 
gas will emit exactly as much radiation to the sensor as 
it absorbs from the blackbody, and the sensor will see the 
same radiation it would see from the blackbody alone. 
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There are other optical methods; however, all of these are 
subject to a "long path length" error. The path length 
includes two boundary layers at reduced temperatures , with 
temperature gradients as great as 6000 P°/inch, and a mass 
of "core" gas which may or may not be homogeneous. How the 
"long path length" instrument weights these different den- 
sity and temperature regions is not known, and this adds 
uncertainties to data interpretation. Some pyrometers 
approximate the mass average, and some the "free stream" 
temperature, and in the absence of an understanding of the 
flow processes and heat transfer, a rational decision is 
difficult. Ghandhi (48) gives an interesting discussion 
and summary of these matters. 

A truly "local" gas temperature measurement was needed for 
this study, and under restrictions in frequency and gas 
behavior, the small-wire transducer can approach this stan- 
dard and give at least a good qualitative, if not a quan- 
titative, understanding of the gas temperature behavior. 

SMALL WIRE INSTRUMENTS 


There are four general categories of instruments employed 
in the study of fluid properties which use a small diameter 
wire as the transducer: the thermocouple, the constant 

temperature anemometer, the constant current anemometer, 
and the resistance thermometer. The principal differences 
between these four categories are the magnitude of the wire 
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current and the philosophy applied to the output signal. 

It is instructive to consider all four categories briefly 
to understand the heat transfer principles involved. 

SMALL WIRE THERMOCOUPLE 

The small wire thermocouple can be used to measure tran- 
sient gas temperatures if the fluctuation frequency is 
relatively low. A thermocouple constructed of 0.001 inch 
diameter wire will have about a 7 millisecond response 
time (to 63.2$ of a step change) in a high velocity gas 
flow. Thermocouples of this sort are available from Heat 
Technology Laboratories*. Omega Engineering, Inc.* makes 
a 0.0005 inch thermocouple; thermocouples as small as 
0.001 inch can be built using resistance welding or spot 
welding joining techniques. However, the frequency re- 
sponse of these instruments limits their applicability to 
speeds less than 100 or 150 rpm (about 2 Hz); they would 
clearly be unsuitable for this investigation. 

CONSTANT TEMPERATURE AND CONSTANT CURRENT ANEMOMETERS 

The constant current anemometer uses a resistance wire in 
a bridge circuit with a high impedance dc power supply to 
insure a constant probe current during small changes in 
probe resistance. The wire is heated with a dc current of 


*See List of Manufacturers and Addresses 
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from 30 to 300 milliamps; when the heated wire is 
cooled by a fluid stream, its electrical resistance de- 
creases with its temperature. The energy dissipation and 
fluid velocity are related by the famous law derived and 
verified by L. V. King (49) of McGill University in 1914 : 

I 2 R = (R - R q ) (A + B /Re) (4) 

where A and B are functions of the gas properties and wire 
dimensions. Re is the Reynolds Number associated with the • 
wire, I the current, R the hot wire resistance, and R the 
room temperature wire resistance. King’s Law was subse- 
quently found to be correct for wide ranges in Re and 
accurate for ordinary wires in the velocity range 0 to 200 
ft/sec. The use of this instrument for time-varying velo- 
city requires electronic compensation and essentially con- 
stant fluid temperature. 

The constant temperature anemometer uses the same basic 
circuitry as the constant current device, except that the 
bridge unbalance signal is passed to a control amplifier 
which automatically adjusts the heating current in the wire 
(still on the order of hundreds of milliamps) to restore 
the bridge to balance. Thus, as long as the amplifier is 
able to bring the bridge near balance, the resistance is . 
maintained constant, and the current in the hot wire is a 
direct measure of the stream velocity at any time (when in- 
terpreted by King's Law). 
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RESISTANCE THERMOMETER 

The resistance thermometer also consists of a fine wire 
immersed in a gas (although the wire must be finer and the 
gas cleaner than in the previous devices), a bridge, a 
power supply (constant current source), a calibration meth- 
odology, d.c, amplifiers, and a display device. As the 
temperature of the gaseous medium changes, the. fine wire 
tries to remain in thermal equilibrium with the gas; as 
its temperature changes, its resistance changes in a well 
known fashion. The resistance changes cause a change in 
bridge unbalance (the bridge never need be in exact bal- 
ance), and this voltage is amplified, displayed, and re- 
corded. Since the wire is very small, its thermal inertia 
is small. It is insensitive to gas velocity except insofar 
as greater gas velocities increase the convective heat 
transfer to the wire, enabling it to remain in closer ther- 
mal equilibrium with the gas. The electrical current 
through the wire must be so small as to make the energy 
dissipation through I 2 R losses negligible. It is in this 
respect, the power dissipation in the wire being negligible 
and the device being insensitive to King’s Law, that the 
resistance thermometer differs from the hot wire anemometer. 
The wire is not "hot" as viewed by its environment. 

Circuitry 

A simplified circuit diagram for the resistance thermometer 


172 


is presented in Pig. 7. The bridge, whose unbalance is 
displayed and recorded, is an equal-arm constant-current 
bridge. 

As the input impedance of the amplifiers and oscilloscope 
is on the order of 100 megohms, negligible current flows 
in the output circuit. The nominal wire resistance and 
dummy arm resistance are about 10 ohms, while .the upper 
arms are 47,000 ohms each; all resistors are 1# low-noise 
film resistors, and all pots are wirewound. The dummy and 
transducer arm currents are fixed by the 47K arms and the 
bridge voltage, and are held constant at a maximum of 0.491 
ma (compared to hundreds of ma for an anemometer) . Since 
wire resistance varies linearly with temperature in the 
range 0°P, to 800°P, it is apparent that a change in bridge 
unbalance is very closely proportional to a change in wire 
temperature. With proper scaling, calibration, and ampli- 
fication, the displayed signal may be taken for wire tem- 
perature. If the wire is small enough for the highest im- 
portant harmonic of the gas temperature fluctuation, the 
displayed signal can be taken as gas temperature. 

During calibration and set-up, switching transients in the 
probe current were thought responsible for a rash of wire 
failures. A ten-turn 300K ohm rheostat was installed in 
series with the power supply and bridge; bridge voltage 
could be continuously varied from 1.8 to 22.4 volts, and 
all subsequent switching was done at the lowest voltage. 



Fig. 7 Resistance Thermometer Circuitry 
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All data were taken with full 22 A volt bridge voltage. 
Problems with switching transients were also noted by 
small-wire experimenters in England, and essentially the 
same solution was employed. 

Fabrication 

Small diameter resistance wires have been used, by experi- 
menters in "homemade" resistance thermometers for the study 
of low frequency transients in gas temperature for over 90 
years. The major drawback to using this type of instrument 
has been its fragility, since extremely small wires must be 
used to provide a high enough frequency response. In very 
recent years, very small wire hot wire anemometer probes 
have been available commercially, but they are generally 
designed for use in the study of fluctuating fluid velo- 
cities at constant temperature. No other experimenters 
are known to have applied these probes to the study of 
fluctuating gas temperatures. 

Since anemometer probes are available in wire diameters in 
many steps from five microns (one micron = 10 ^ meters) to 
less than one micron, in tungsten, platinum, or platinum/ 

10/5 rhodium, and in various probe styles, they seem to be 
natural transducers for the study of changing gas tempera- 
tures at a point. However, there are a number of distinct 
disadvantages and limitations to their use in this applica- 


tion: 
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1) The environment must be clean and free from 
oil droplets, dust, or condensed phases. 

2) Local velocities may not exceed a few hundred 
ft/sec. (see Appendix B) , the mean temperature 
should not exceed about 300°P, and peak temper- 
atures should not exceed about 800°F. 

3) Initial cost is in excess of $50 per probe, 
with a $20 average for wire replacement. 

4) It is difficult to seal the probe bodies 
against compression pressures above 200 psi. 

5) Wire length is short and electrical resistance 
is very low, making the design of a satisfac- 
tory bridge difficult. 

To counter these difficulties, it is sometimes necessary 
to design a custom resistance wire mounting to the system 
specifications . 


The custom wire mountings initially used in this project 
were constructed on a half-inch thick, one-inch diameter 
cast iron instrument plug with a 1-14 UNF thread. Fig. 8 
shows two views of such a plug with a centrally located 2-56 
UNC threaded hole for a surface thermocouple. The resis- 
tance wire used is Wollaston process wire; this wire is 
supplied as a one mil (25.4 micron) silver wire with a one 
to five micron core wire of platinum/10# rhodium. The 
silver sheath is etched away after mounting, leaving the 
precious metal core as the active resistive element. Two • 

20 gage constantan wires were used with a 3/8 to 5/8 inch 
spacing as wire supports. Constantan was used because it 
has low thermal conductivity and heat capacity while re- 
taining high electrical conductivity. 
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A OLD MOUNTING METHO 




B NEW MOUNTING METHOD 


CONSTANTAN SUPPORT 



Fig . 8 hand-made instrument plugs 
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The most difficult problem in probe construction is making 
a rigid, electrically conductive, low heat capacity connec- 
tion of resistance wire to supports. If the resistance 
wire diameter is less than .0005 inches (12 microns), it 
is practically impossible to solder, spotweld, or otherwise 
effect a satisfactory metallic bond of wire to constantan 
support because the wire is easily melted. Cold solder 
joints tend to soften under high compression temperatures, 
and flux flows along and coats the wire. Silver solder 
joints were attempted using the smallest torch commercially 
available*, with a butane -oxygen flame 1/2 millimeter by 
1 millimeter; the wires melt or deform at the conditions 
required to melt the solder. In recent years, micromani- 
pulative and microwelding equipment capable of handling 
wires as small as one micron has become available. This 
equipment is used by commercial suppliers of hot wire ane- 
mometer probes; it is, however, very expensive • and requires 
a talented and patient operator. 

Silver-impregnated epoxy resin* (designed for use in print- 
ed circuit board repair) was used for a time by the author 
to hold the wires to the supports, but the epoxy resisti- 
vity is quite high, and the epoxy joint constituted a very 
high resistance (about 400 ohms) in series with the active 
length of wire (3 or 4 ohms). As this series resistance 
was temperature-sensitive and the silver impregnation des- 


*See List of Manufacturers and Addresses 
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troyed most of the bonding power of the epoxy, the re- 
sistance thermometer performance was erratic; this bond- 
ing method was therefore discarded. 

The best solution for the mounting problem was later 
corroborated by a private communication with G. J. Schott, 
associate of B. H. Schultz*, dean of resistance thermometer 
experimenters. Fig. 8B shows the construction of an in- 
strument plug incorporating the solution. The wire is 
stretched and suspended, with its silver sheath intact, 
over the plug and its constantan support posts. The 
assembly is held on a special jig on the stage of a Leitz 
stereoscopic microscope. The wire is lowered until it i.s 
stretched across the centers of the constantan post tops 
and is placed under moderate tension. This is done be- 
cause the core wire is kinked slightly inside the silver 
sheath. Using fine dental tools and a 50x gain on the 
microscope, insulating household two-component epoxy is 
daubed onto the wire-support junctions and allowed to dry 
and harden for 24 hours. The wire excess outside of the 
active length between the support posts is wrapped around 
two 2-56 brass machine screws which serve as electrical 


*Schultz was a brilliant physicist and experimentalist 
employed by Phillips in Holland. According to G. J. Schott 
Schultz straightened his wires and reduced their diameters 
by rolling them between glass slides separated at one end 
by a gramophone needle. The wires were used in Stirling 
cycle devices using air and light gases as working fluids. 
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terminals. Thus the constantan posts serve only as 
mechanical supports and to define the active length, and 
the machine screws carry the electrical signal. Wire slack 
is allowed between support post and screw, and the screw 
is tightened down and sealed with nonconductive epoxy. 

After this seal dries, the etching operation can commence. 

As mentioned earlier, the wire is mounted in the form of 
a small platinum/10# rhodium core inside a silver sheath. . 

The silver is etched away with a drop of dilute nitric 
acid: a 40# U^O / 60# HNO^ mixture worked best and pro- 
vided an etching time of 25 to 35 seconds. During the 
entire mounting and etching operation, the central surface 
thermocouple was replaced with a steel 2-56 machine screw. 

As the etching mixture also attacks constantan, cast iron, 
and epoxy, it is necessary to protect these surfaces with 
a non-etchable coating. Schultz used nail polish. The 
author used paraffin, which was applied with a small, heated, 
hollow brass tube to the support posts, epoxy joint, and in 
a circular "wall" which defined and concentrated the etch- 
ing droplet. One or two drops of etch mixture are applied 
with an eye-dropper; the etching action is watched through 
the stereo microscope. A skin of paraffin forms on the 
etch droplet and will grip the wire and tear it from the 
supports when the etchant is drained away. This skin is 
dissolved by adding 1/2 drop of xylene to the nitric acid 
mixture at the end of the etching period. The fluid is 
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then quickly removed with a 25 gage hypodermic needle. 

The surface is flushed with two or three drops of distilled 
water and the protective paraffin is dissolved with drops 
of xylene. Residual paraffin is removed with curved dental 
tools; this is a very delicate step. 

Using this technique, over a dozen probes were constructed 
with wires from 1/8 inch to 5/8 inch long and penetration 
depths of 1 to 5 millimeters. Several multiple mounts 
were made with wires parallel or crossing at different 
penetration depths. One plug was constructed in which the 
penetration depth could be continuously varied from near 
zero to 5 millimeters. 

Platinum Plating Technique 

A possible alternative to the fragile small wire was a 
platinum-plated quartz wire. According to a monumental 
work by H. H. Lowell and N. Patton (50), under certain 
conditions of frequency and geometry a solid cylinder of 
insulating material with a thin metallic coating will have 
a frequency response superior to a solid metal wire of the 
same size. 

Quartz fibers of .0004 inch (9 micron) diameter were ob- 
tained from the General Electric Company*. These filaments 
have a tensile strength in excess of 500,000 psi and could 


*See List of Manufacturers and Addresses 
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be mounted by silver soldering at 1100°C. 

A plating solution was obtained* which contained organic 
gold and platinum compounds in a solution of volatile oils. 
When fired at 1022°P, the solution will coat a fuzed quartz 
fiber with a 91 % pure conductive film 0.000005 inches (0.1 
micron) thick with a resistivity of 35 micro-ohms/centimeter. 

Preparations were underway to attempt the difficult mount- 
ing and firing operations when it was discovered that a 
manufacturer of anemometer probes had the capability of 
constructing probes in which the wire was replaced by a 
25.4 micron hollow glass tube with the fired precious metal 
coating 0.5 microns thick. As this size tube would signi- 
ficantly disturb the flow patterns and hold a large boun- 
dary layer, inquiry was made whether the firm** could con- 
struct a probe with a 10 micron or smaller tube. After 
lengthy consultation, it was determined that the number 
of fired coatings required to produce a reasonable thick- 
ness of metal increased inversely as the square of the dia- 
meter and became unreasonable at 20 microns and smaller. 

The basic principle of this type of system is that at high 
enough frequencies the core becomes thermally dissociated 

from the coating and remains at constant temperature; the 

> 

*See List of Manufacturers and Addresses 

**ThermoSy stems, Inc.: see List of Manufacturers and 

Addresses 
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thin conductive coating follows the environment tem- 
perature. At lower frequencies the entire laminated cylin- 
der responds as a solid wire (i.e., very slowly). The 
Lowell reference (50) previously mentioned indicated that 
for the geometries commercially available, the core ma- 
terial would not dissociate itself thermally below 100 cps 
(6000 rpm); for these reasons the coated quartz fiber 
probe was not used. 

Phase Lag and Attenuation 

This discussion will restrict itself to very small wires 
(.001 inch and smaller) under the condition of negligible 
self-heating. When a small wire is exposed to a gas whose 
temperature is changing, the thermal inertia. of the wire 
and energy absorption and liberation in the boundary layer 
surrounding the wire will prevent it from exactly following 
the temperature change. The wire will tend to lag in time 
(a phase lag) and experience a smaller temperature swing 
than that of the driving potential (an amplitude attenua- 
tion) . The smaller the wire for a fixed environmental 
history, the smaller the boundary layer around it, the 
lower the inertia, and the smaller the phase lag and am- 
plitude attenuation. It was the intention of the author 
to use smaller size and faster response instruments than 
previous investigators, since the methods of compensating 
for lag and/or attenuation are on the whole not very satis- 
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factory. 

An interesting comparison of resistance thermometer and 
infrared detector (see page 39 for description of the 
latter) was presented by H. J. Wienke (51) of the Univer- 
sity of Wisconsin. Wienke compared infrared null method 
gas temperature data of K. C. Tsao and co-workers (52), 
taken in an unfired CPR diesel engine, with his data, also 
taken in a CFR head and crankcase with tungsten wire 
sizes from .001 inch to .00015 inch. Wienke extrapolated 
peak cycle temperature ana peak crankangle linearly with 
wire diameter, taking the infrared null data to represent 
a zero diameter wire. This linear extrapolation suggested 
that the finest size resistance wire tested attenuated the 
peak gas temperature by about 150°F and lagged the peak 
crankangle by some five or six crankangle degrees. 

Several objections may be presented with regard to the 
Wienke data and its interpretation. Objections to the 
infrared null method on the basis of long path length 
errors are noted elsewhere in this document . The signal 
is integrated along the entire path length, which includes 
two cooler boundary layers. Comparison of the "weighted" 
gas temperature with a locally measured gas temperature 
must be made with caution. 

In the second place, all the data were taken with tungsten 
wires in an unmodified CFR compression-ignition engine with 
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oil-lubricated cylinder walls and crankcase. Experience 
shows that within a few dozen cycles from starting motor- 
ing, the air contains minute oil droplets and the metal 
surfaces of the compression chamber are covered with an 
oil film. This oil accumulation was the primary cause for 
the very small wire breakage noted by the experimenter; 
Schultz (^0) had similar troubles during experiments in a 
Stirling cycle device. The effects of oil and dirt accum- 
ulation on the wire are noted by H. P. Grant (53) in an 
American Society of Mechanical Engineers presentation on 
anemometry : 

"Foreign matter in the flow gradually accumulates on 
a hot wire and ultimately changes its calibration . 
(and response) by altering the heat transfer rate. 
...At high (fluid) velocity the heat transfer rate 
is always reduced. ...Without air filtering, clean- 
ing will be required every few hours at velocities 
below 100 ft/sec., but may be required every few 
minutes at a velocity about ^00 ft/sec. to avoid 
large corrections for dirt (and oil) accumulation. ' 
Tungsten wires must be cleaned with . . . K^Cr-O,. sat- 
urated with H 2 S0^. !! (parentheses added ) d ' d 1 

With reduced heat transfer rate the dynamic response of 
the wire is reduced: lag and amplitude attenuation are 

increased. Thus, wires in a dirty or oily atmosphere 
will experience response difficulty. 


A final consideration is that extrapolations of peak tem- 
perature and crankangle should not be made linearly with 
wire diameter. The classical wire time constant varies 
with the square of the diameter; the variation of mean heat 
transfer coefficient about the wire with diameter must also 


185 


be taken into account. Based on data presented by 
Schultz (40) and conversations with several professional 
persons knowledgeable in small wire work, the index 1.7 
was chosen as more appropriate. When Wienke's original 
data is extrapolated with respect to d^*^, a phase lag 
close to one crankangle degree and amplitude diminutions 
of 20 F° (at 600 rpm) and 30 F° (at 1200 rpm) are suggest- 
ed for the 3.81 micron tungsten wire. This is' much closer 
to analytical predictions of lag. 

If a wire is in equilibrium with an environment at tem- 
perature T q and the environment suddenly steps to T-, , then 
if we ignore capacitance in the boundary layer about the 
wire and assume the film coefficient to be a constant, the 
wire temperature at any time is: 

T = T x - (T x - T q ) exp (-hAt/Wc) ( 5 ) 

and the time constant of the wire is: 

6 = Wc/hA (g ) 

where W is the weight of the wire, 

c is the specific heat of the wire material 
h is the film coefficient, and * 

A is the wire surface area. 

If we take h as the wire diameter— dependent value calcu- 
lable from Schultz' (40) data, and define the phase lag 
as the time period for 63 . 2 % response of the wire, then • 

9 gives the phase lag of a platinum wire. These 
calculations ignore capacitive storage and pressure work 
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in the wire boundary layer, and should be taken as an 
approximate guide only. The analogous plot for tungsten 
is little different from Pig. 9. 

Self-Heating of Wire 

It is fundamental in the operation of constant current 
and constant temperature anemometers that the ■wire dissi- 
pate relatively large I 2 R losses to its environment in 
order that King's Law might be used to calculate the local 
fluid velocity. The ideal resistance thermometer, however, 
would exchange no energy with the gas (and consequently 
the difference [T wire - T gas ] would be zero). Since a 
small current must pass through the wire for sensing pur- 
poses, it is necessary to set an upper current limit at 
which the self-heating is no longer insignificant . Typi- 
cal operating conditions with a D.I.S.A. 55A35 < subminia- 
ture probe are as follows: 

Diameter of wire = 3.81 microns 
Length of wire = 2 millimeters 
Resistance of wire = 3 ohms 
Mean film coefficient ■ 30 Btu/hr ,°P.ft 2 

The last figure represents a conservative estimate of the 
mean film coefficient around a very small cylinder under 
compression or expansion gas velocities. The wire area 
available to dissipate heat is: 

Area ®fiDL *»7t(3.8l microns) (2mm) = 2.58 x 10“^ft 2 . 
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If we allow a maximum wire-to-gas temperature differ- 
ential of 1 F°, the maximum allowable energy flux is: 

Power = hA(T w -T g ) = (30)(2.58 x 10" 7 ) (1)/(3.413 Btu/watthr) 
Power » 2.26 microwatts 


This means: 



7 2.26 microwatts 
3 ohms 


0 . 87 ma 


(7) 


Since our wire current was limited to .48 ma, we have 
operated within the allowed limit. Under the specified 
conditions, a current of 0.48 ma would produce a tempera- 
ture rise of 0.3 P°, a negligible amount compared to a 
600 F° gas temperature swing. 


Strain Gage Effects 


So-called strain gage effects could in principle occur 
when high fluid velocity causes the wire to bow between 
its supports. Stretching of the wire would cause its 
resistance to change; this would look to the electronics 
like a fluid temperature change . 


Weske (54) discusses this effect and suggests that in 
past experience it has caused little trouble. The fun- 
damental strain gage relation is: 




GF 


( 8 ) 
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Weske uses a Gage Factor of 1,7 for tungsten and 2.15 
for cupronickel; we will take a general GF of 2.5 as a 
conservative estimate. Deflections observed under the 
microscope while the probe is exposed to a high velocity 
air stream suggest a AL/Lq of at most two or three per- 
cent. With formula (8) above, 

AR/R q = C2.5X.03) = .075 

In a typical data run, the wire temperature swing is on 
the order of 600 F° . With a temperature coefficient of 
resistance of .0039 ohms /ohm C° for platinum, the wire 
resistance swing during a data cycle would be calculated 
as. follows: 

R(T) = R q + R q o(T - T q ) 

AR = R(T) - R q = RqOiAT 
AR 

— = cxAT = (. 0039/0°) (600F°)(5C°/9F°) 

0 

^- 1.60 (9 ) 

Thus the predicated strain gage effect would be less than 
5/» of the total resistance swing. Experimentally, when the 
probes were exposed to a jet of 70 psi air exhausting 
through a .07 inch diameter hole, no perturbations in the 
bridge output signal could be observed. 

According to Weske, a resonant vibration in the wire is 
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typically started and maintained by the flow, with 
accompanying cyclic stresses. The frequency of this wire 
vibration has been experimentally determined to be: 

f - 0.22 V n /d (10) 

where V n is the normal component of stream velocity, 
and d is the wire diameter in feet. 

I 

(■ 

For an engine with V n on the order of 100 'ft/sec and with I 

3.81 micron wire (the size most frequently used), the 

i 

frequency of the vibration, and of the cyclic stress, is | 

1.8 megaHz. The frequency response of the D.C. ampli- 
fiers is such that this signal will be greatly attenuated | 

and will cause no trouble. \ 

t 

• I 

f 

Other Considerations | 

Three other considerations periodically mentioned in the 
literature pertaining to small resistance wires are aging, 

wire support end effects, and slip flow conditions. j 

i 

! 

I 

Aging of the wire means changes in its room temperature \ 

resistance due to grain growth, grain boundary changes, j 

and inelastic stretching during the high-temperature por- 

I 

tion of its operation. This effect is relatively long- l 

term and is compensated for by recalibration of the wire 
prior to each data run. 


Wire support end effects refer to the influence of the 
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wire supports as heat sinks on the surrounding fluid. 
Aftalion (55) claims that fine wire butt-welded thermo- 
couples are five times as sensitive to end conduction 
effects as are resistance thermometer wires. This effect 
is deemed negligible for wire L/d ratios in excess of 
certain somewhat arbitrary limits. Hunt (56), as reported 
by Meyer (57) > calculates a conduction error of less than 
12 with an L/d of 52 and "heavy lead wire" supports. 

Grant (53) gives a lower limit of 100 for small, pointed 
supports; Thermo-Systems, Inc. (58) states that an L/d of 
at least 400 is common. For larger support sizes, Schultz 
(40) maintains an L/d of 2000. An outline of these con- 
siderations is given by references 59 and 60. For two 
styles of Wollaston wire probes, the author fixed the 
ratio at 1900 and 3200; for the D.I.S.A. 55A35 anemometer 
probes it was 524, removing support end effects from con- 
sideration. 

One must consider the possibility of slip flow occurring 
for very small wires, since a transducer in the slip flow 
or transition regions will not see its environment as a 
continuum. According to Otis (61), the transition, slip 
flow, and "continuum" regions could be approximated by 
the limits: 

transition flow: 0.1 < - <10 

10 < £ < 100 


slip flow region: 


(11 A,B,C) 


continuum region: 100 < — , 

where d is the significant transducer dimension (i.e., 
wire diameter) and X is' the gas molecule's mean free path 
A convenient approximate formula for X is: 


X (in centimeters) 


5 

P (microns Hg) 


( 12 ) 


Then over a range of 20 psi to 400 psi the gas mean free 
path varies between 0.048 microns and 0.002 microns. The 
d/X ratio for the 3*81 micron wire used varies from 80 
to 1600. Thus on the whole the flow past the wires used 
was continuous; however, for smaller wires and/or lower 
pressures, care would have to be taken in interpreting 
the data. 


MULTIPLE WIRE METHODS 

In an attempt to circumvent the problems of phase lag 
and amplitude errors caused by wire and boundary layer 
inertia, some experimenters have proposed multiple wire 
schemes. in which simultaneous temperature measurement 
with two or more transducers of differing inertia would 
yield a solution for instantaneous gas temperature. 

Pfriem (62) made use of energy considerations in an early 
work to correct the temperature indications of several 
sizes of resistance thermometer. Rounthwaite (63) em- 


ployed a two-wire resistance thermometer in the study 
of temperatures in explosions of gases. 
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Aftalion (55) of Sulzer used two very large (50 and 100 
micron) butt-welded thermocouples to study exhaust gas 
temperatures in a huge stationary diesel engine. Recor- 
ding his data on a Siemens Oscillograph (the frequency 
response of that instrument must have restrained the 
useable cyclic frequency of the engine to a few hundred 
rpm) and correcting for radiation, convection, and capa- 
citive inertia by means of a nomograph, Aftalion came up 
with some unusual calculated gas temperatures. A sample 
of his calculated and measured temperatures is reproduced 
at right. 


According 
to Afta- 
lion, be- 
cause of 
the com- 
bustion, 

"a fine 
deposit of 
soot is 
expected 



(on the wires) after a very short period of operation". 
This soot will increase the local heat transfer coeffi- 


cient and the area available for convection (see page 56 
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for the effect of dirt and oil on calibration) j this 
effect was not considered by Aftalion. 

Aftalion measures wire temperature swings of 2^5 P° and 
380 F° with the two large wires and from the slopes of 
these data computes a gas temperature curve with a swing 
of 930 P° . An accuracy figure of 2 % is claimed, and various 
physical interpretations are given the oscillations in the 
computed curve. 

A simplified two-wire technique is advanced by W. E. Meyer 
(57) in the Society of Automotive Engineers Transactions. 
Equating the convective heat transfer to a thermocouple 
wire to its rate of change of internal energy, for each of 
two thermocouples, A and B (subscripts), 

(pcd/4)(T* - T) - h ( T - 1/2[T' + T])Ax (13 A,B) 

gaS 

where p is wire density, 

c is wire specific heat, 
d is wire diameter, 

h is the film coefficient of heat transfer, 

T _ is the average gas temperature for the 

gcio 

period At, and 

T, T' are junction temperatures at beginning 
and end of the time period At. 

This equation is reproduced here because the many print- 
ing errors in the Transaction version may cause confusion. 
Designating the two thermocouples by the subscripts A and 
B, we write Nusselt's expression for each thermocouple: 
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( 14 ) 


(15) 


Taking the ratio of these relations to eliminate common, 
but unknown gas properties, we get: 




1-n 


( 16 ) 


Combining Eqn. 16 with Eqns. 13 A, B, and solving for gas.' 
temperature, we get an expression for gas temperature as 
a function of the unknown n and a number of known quan- 
tities: 


Vs - f(n) - 


(17) 


Meyer proposes to solve for the unknown T from among 
the four unknowns T gas , h A , h B , and n with the simultan- 
eous solution of three independent equations (note that 
[17] above is not independent): 


h A * f l<Vs> 


hr, * f~(T ) 

B 2 gas 


(or Eqn. 13 A), 
(or Eqn. 13B), 


n > f 3 (h A> h B ) 


(or Eqn . 16) . 
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The technique given is to assume an n (aided by values 
from McAdams' book [64], which are based on early data 
of L. V. King, the formulator of King's Law of Anemom- 
etry), calculate a T gas , solve for h A and h Q with Eqns. 

13 A, B, and verify the assumed n with Eqn. 16. This 
calculation is clearly not possible; any choice of n at 
all will work equally well. Even then, it is not clear 
that the Nusselt relation, measured in steady-state gas 
flow, is valid for cylinders only a hundred times a mean 
free path in diameter in a very unsteady gas environment. 


Somehow using this method, Meyer calculates a smooth gas 
temperature from data taken in a motored diesel engine at 
100 and 150 rpm and finds it to be only a few degrees above 
the indicated thermocouple temperatures at top dead center 
for 20 and 13 micron thermocouple wires. At these fre- 
quencies (less than 3 Hz) the thermocouple errors are so 
small (3# at most) that the correction, regardless of how 
it was obtained, is trivial, and does not affect the vali- 
dity of the conclusions Meyer drew from the data. 

A three-wire method expounded by Benson is being used by 
G. A. Morgan (65) of the University of Adelaide in Aus- 
tralia in a work similar to this one. Benson's differen- 
tial equation reduces to: 


T 

gas 


T 

wire 


+ d 


2-m 

dT wire] 

[w m . 

Pw C 


l dT j 


4AXg 

m 


( 18 ) 
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with the various constants being evaluated from the 
three-wire data. Benson comments that d^” m is not sensi- 
tive to small changes in m, as m is much smaller than 2. 

An extension of Benson’s work by Daft shows that three 
wires enable m and all of the bracketed terms to be treated 
as unknowns and eliminated in the solution for T 

gas 

Multiple wire methodologies generally depend on the slopes 
of experimental data and are subject to scatter; in addi- 
tion, they ignore the capacitive and pressure work pheno- 
mena thought to exist in large-wire boundary layers and 
their effect on wire signal lag and attenuation. 

In view of these uncertainties, multiple wire techniques 
were not used in this work, and emphasis was placed on 
obtaining repeatable signals from the finest diameter re- 
sistance wire possible, consistent with the requirement of 
moderate durability. 


\ 


VI . CALIBRATION 
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PRESSURE RECORDING SYSTEM 

The pressure pickup used was a Kistler 601H piezoelectric 
transducer. A Li-Draper catenary diaphragm (strain gage) 
pickup was originally used, but the Kistler diameter is 
less than that of the Li, and using it allowed more flex- 
ibility in instrument placement. 

The Kistler pressure pickup was calibrated under static 
conditions in a dead weight tester and its output was 
found to be essentially linear over the range —11 psig to 
600 psig; with an average output coefficient of 1.078 pc b/ 
psi over this range, the transducer output was about 8 % 
above the nameplate output of 1.000 pcb/psi. The pickup 
output was calibrated on the display by noting the zero 
gage pressure level on the oscilloscope before- and after 
the data run with the charge amplifier in short time con- 
stant operation. 

SURFACE THERMOCOUPLE 

The surface thermocouple was calibrated as follows: placed 

in & lucite block which protected all rustable surfaces, 
the thermocouple was connected in series with a cold junc- 
^ion and a Leeds and Northrup potentiometer and was immersed 
in stirred distilled water of different temperatures from 
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0°P to 212°P. Fig. 10 shows the result of this calibra- 
tion; the outputs were closely linear, and the most repre- 
sentative slope was taken to be .0298 mv/P° for iron- 
constantan thermocouples and .0170 mv/P° for iron-nickel 
thermocouples. Eleven thermocouples were calibrated, in- 
cluding four pre-1960 and four newly purchased iron-nickel 
thermocouples. 

Knowing the thermocouple gain, d.c. amplifier gain, and 
Tektronix oscilloscope gain (individually and collectively 
calibrated), it was necessary to know the temperature at 
one time on the oscilloscope’ display or on an exposed drum 
camera film to fix all the head surface temperatures. 
Previous investigators used automotive breaker point assem- 
blies connected to camshaft or other geared-down power 
takeoffs. In keeping with this philosophy, the author 
spent some time designing and constructing a similar assem- 
bly to run at the crankshaft speed and provide adjustable 
intervals in each engine cycle during which known voltages 
could be superimposed upon the various traces . Because 
this mechanism had to run up to 3600 rpm (as opposed to 
800 or 900 rpm for previous CFR camshaft-driven setups), 
and because the present signal levels were very low, find- 
ing the signal amidst the noise was difficult. 

A simpler technique suggested by G. Ebersole (66) was then 
adopted, wherein a Leeds and Northrup type 7787^6 potentio- 


SURFACE THERMOCOUPLE OUTPUT 
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meter was used to measure the time-averaged head surface 
temperature; this temperature was used as the first term 
in the Fourier breakdown of the surface temperature trace. 
This technique is much easier and faster to use than the 
previous method, and it doesn't require sending low level, 
high frequency signals through automotive breaker points, 
significantly improving the signal/noise ratio. 

RESISTANCE THERMOMETER 

The resistance thermometer calibration posed some problems. 
The probe was sealed inside the compression chamber with 
piston at bottom dead center and hot air was circulated 
about the cylinder and head from an external blower and 
bank of heaters, the air system of which is shown of Fig. 1. 
When the system had reached equilibrium at a known temper- 
ature, the bridge output for a given bridge voltage was 
recorded. Many times, however, an expensive resistance 
wire which could withstand compression and expansion heat- 
ing and gas velocities would break under this gentle ex- 
ternal heating. Greater success was achieved when the 
heating method was changed: hot air was blown through the 

makeup holes in the lower cylinder wall, across the wire, 
and exhausted through a port in the head. It is postulated 
that uneven wire support expansion snapped the wires in the 
previous method. 


Fig. 11 shows bridge output as a function of wire tempera- 


bridge output calibration / 



Fig. 11 


Resistance wire Calibrati 


ion 





ture for a typical D.I.S.A. probe; the relationship 
was linear through the 70°F - 190°F test range and was 
assumed to be linear over the whole range of operation. 

Resistance thermometer output was calibrated on the dis- 
play by noting the vertical location of the room tempera 
ture line on the oscilloscope just before operation and 
setting this abscissa equal to room temperature. 


